Quantum Impurity Problems in Condensed Matter Physics [□] 

Impurities are ubiquitous in condensed matter. Boundary Conformal Field Theory (BCFT) pro- 
vides a powerful method to study a localized quantum impurity interacting with a gapless continuum 
of excitations. The results can also be implied to nanoscopic devices like quantum dots. In these 
lecture notes, I review this field, including the following topics: 

I. General Renormalization Group (RG) framework for quantum impurity problems: example of 
simplest Kondo model 

II. Multi-channel Kondo model 

III. Quantum Dots: experimental realizations of one and two channel Kondo models 
. . , IV. Impurities in Luttinger liquids: point contact in a quantum wire 

0^ ' V. Quantum impurity entanglement entropy 

' VI. Y-junctions of Luttinger liquids 

VII. Boundary condition changing operators and the X-ray edge problem 

(N . 
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CD I. QUANTUM IMPURITY PROBLEMS AND THE RENORMALIZATION GROUP 

Ql .... n 

, A remarkable property of nature, that has intriqued physicists for many years, is universality at critical points, [ij 
An impressive example is the critical point of water. By adjusting the temperature and pressure, a critical point is 
[ reached where the correlation length diverges and the long distance physics becomes the same as that of the Ising 

1 I '■ model. A microscopic description of water is very complicated and bears very little connection with the Ising model; in 

I particular, there is no lattice, no spin operators and not even any Z2 symmetry. Nonetheless, various experimentally 
• measured critical exponents appear to be exactly the same as those of the Ising model. Furthermore, the best 
' description of this universal long distance behaviour is probably provided by the ip^ field theory at its critical point. 

1 t , Our understanding of universality is based upon the RG. For a system at or near a critical point with a diverging 
d correlation length, it is convenient to consider an effective free energy (or Hamiltonian) , used only to describe long 
^ , distance properties, which is obtained by integrating out short distance degrees of freedom. It is found that the same 

long-distance Hamiltonian, characterizing an RG fixed point, is obtained from many different microscopic models. 
These fixed point Hamiltonians are universal attractors for all microscopic models, [l] 

The same features hold for quantum models of many body systems at low temperature. In many cases such models 
^ ' can exhibit infinite correlation lengths and vanishing excitation energy gaps. In this situation one again expects 
universality. Important examples are the Fermi liquid fixed point for interacting electrons in 3 dimensions (D=3),[2| 
its cousin, the Luttinger liquid fixed point in D=l and various models of interacting quantum spins.[3, 13| 

In these lectures, I will be concerned with a single quantum impurity embedded in such a critical system. The 
' quantum impurity can be of quite a general form, possibly comprising several nearby impurities. If we study its 
[~-«, [ behaviour at long length scales (compared to all microscopic lengths including the spatial extent of the impurity 
> and the range of its interactions with the host) and at low energies compared to all microscopic energy scales, then 
CO universality again emerges. These single impurity models, while simplified, have the attractive feature that such 
powerful methods as BCFT can be used to tackle them. They provide quite non-trivial examples of quantum critical 
■ phenomena and, in some cases, appear to be good descriptions of experimental reality. 
00 . An important example is provided by the simplest version of the Kondo model. [1, 0, This is a model invented 
' to describe a single magnetic impurity (such as an iron atom) in a non- magnetic metal (such as copper). Traditional 
experiments in this field always involve a finite density of impurities, but if this density is low enough, we may consider 
only one of them; technically this gives the first term in a virial expansion in the impurity density. Furthermore, these 
models can be applied to situations where the single impurity is a nanostructure device like a quantum dot. A simple 
^ . Hamiltonian to describe this system can be written: 



o 
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H = J d^fc4"^,„.(fc) + J J -^^f-fi^y, ■ s (1.1) 

Here t/i^^ annihilates an electron of wave- vector k and spin a and is normalized so that: 

{^l;',i^P0}^6pHk-k'). (1.2) 

Repeated spin indices are summed over. e{k) is the dispersion relation for the electrons, which we will usually 
approximate by the free electron form: 



^2 



(1.3) 
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where ep is the Fermi energy. (So this is actually H — fJ.N.) S is an impurity spin operator, of magnitude S. J 
measures the strength of a Heisenberg type exchange interaction between the electron spin density and the impurity 
spin. Usually J > 0. Note that this form of interaction is a (5-function in position space: 



H = d*r 



(1.4) 



Here we have suppressed the spin indices completely. Actually, this model is ultra-violet divergent unless we truncate 
the integral over fc, k' in the interaction term. Such a truncation is assumed here but its details will not be important 
in what follows. The dimensionless measure of the strength of the Kondo interaction is 

A = Jv, (1.5) 

where v is the density of states, per unit energy per unit volume per spin. For free electrons this is: 

rnkp 

v=^^ 1.6 

where kp is the Fermi wave- vector. Typically, A <C 1. 

This model is a considerable simplification of reality. In particular, electrons in metals interact with each other 
via the Coulomb interaction and this is neglected. This can be justified using Fermi liquid ideas. Since A 1, the 
Kondo interaction only affects electrons close to the Fermi surface. The Coulomb interactions become increasingly 
ineffective for these electrons, as can be seen from phase space arguments, after taking into account screening of the 
Coulomb interactions. The free electron Hamiltonian (with an appropriate effective mass) represents the fixed point 
Hamiltonian, valid at low energies. Treating the Kondo interaction as a 5-function is another approximation; a more 
realistic model would give it a finite range. Again, if we are concerned with the long distance, low energy physics, 
we might expect this distinction to be unimportant. The spherical symmetry of the dispersion relation and Kondo 
interaction will considerably simplify our analysis, but again can be seen to be inessential. 

Due to the absence of bulk interactions, the (5- function form of the Kondo interaction and the spherical symmetry 
of e(fc), we may usefully expand the electron operators in spherical harmonics, finding that only the s-wave harmonic 
interacts with the impurity. (See, for example. Appendix A of jsj.) This gives us an effectively one-dimensional 
problem, defined on the half-line, r > 0, with the impurity sitting at the beginning of the line, r = 0. Thus wc write: 

V'g — — ipo{k) + higher harmonics. (1-7) 
V47rfc 



Next we restrict the integral over k in the Hamiltonian to a narrow band around the Fermi wave- vector: 

- A< k - kp < A- (1.8) 

This is justified by the fact that A <C 1. To be more accurate, we should integrate out the Fourier modes further 
from the Fermi surface, renormalizing the Hamiltonian in the process. However, for small A this only generates small 
corrections to H which we simply ignore. Actually, this statement is only true if A is chosen to be small but not too 
small. We want it to be <C A:^. However, if it becomes arbitrarily small, eventually the renormalized A starts to blow 
up, as we discuss below. Thus, we assume that A is chosen judiciously to have an intermediate value. We can then 
approximate the dispersion relation by: 

e{k)fiiVF{k-kF). (1.9) 
We now define the following position space fields: 

.A 

^L/R{r)= / dfce±^">o(feF + fc). (1.10) 

J -A 

Note that these obey the boundary condition: 

Mt,r^O)^^Rit,r = 0). (1.11) 
Furthermore, they obey approximately the anti-commutation relations: 

{VL/fl(x), ^[/^(x')} = 27tS{x - x'). (1.12) 
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[This is only approximately true at long distances. The Dirac (5- function is actually smeared over a distance of order 
1/A. Note also the unconventional normalization in Eq. (|1.12p .] The Hamiltonian can then be written: 



H — —i / dr ( iJJi^-j-ipL — Wr-t'^R I + '"f^Wl{^)t:^l{^) ' ^ + higher harmonics. (1-13) 



The "higher harmonics" terms in H are non-interacting and we will generally ignore them. This is a (1+1) dimensional 
massless Dirac fermion (with 2 "flavours" or spin components) defined on a half-line interacting with the impurity 
spin. The velocity of light is replaced by the Fermi velocity, vp- We shall generally set vp = 1. Note that in a 
space-imaginary time representation, the model is defined on the half plane and there is an interaction with the 
impurity spin along the edge, r — 0. Since the massless Dirac fermion model is a conformal field theory, this is a type 
of conformal field theory (CFT) with a boundary. However, it is a much more complicated boundary than discussed 
in John Cardy's lectures at this summer school. There he considered CFT's with conformally invariant boundary 
conditions. If we set A = then we have such a model since the boundary condition (BC) of Eq. (jl.lip is conformally 
invariant. More precisely, we have a boundary conformal field theory (BCFT) and a decoupled spin, sitting at the 
boundary. However, for A 7^ 0, we do not have merely a BC but a boundary interaction with an impurity degree of 
freedom. Nonetheless, as we shall see, the low energy fixed point Hamiltonian is just a standard BCFT. 

Although the form of the Hamiltonian in Eq. (|1.13p makes the connection with BCFT most explicit, it is often 
convenient to make an "unfolding" transformation. Since ipLit, x) is a function of (t+x) only and ipRi^y ^) is a function 
of {t — x) only, the boundary condition, of Eq. (jl.lip implies that we may think of ipR as the analytic continuation 
of tpL to the negative r axis: 



^Rir)=M-r), {r>0) (1.14) 



and the Hamiltonian can be rewritten: 



dr^l—4>L + VFXi^liO)^MO)-S. (1.15) 

We have reflected the outgoing wave to the negative r-axis. In this representation, the electrons move to the left only, 
interacting with the impurity spin as they pass the origin. 

The phrase "Kondo problem" , as far as I know, refers to the infrared divergent property of perturbation theory, 
in A, discovered by Kondo in the mid-1960's. In the more modern language of the RG, this simply means that the 
renormalized coupling constant, A(i?), increases as the characteristic energy scale, E, is lowered. The "problem" is 
how to understand the low energy behaviour given this failure of perturbation theory, a failure which occurs despite the 
fact that the original coupling constant A <C 1. The /3- function may be calculated using Feynman diagram methods; 
the first few diagrams are shown in Fig. ([T]). The dotted line represents the impurity spin. The simplest way to deal 
with it is to use time-ordered real-time perturbation theory and to explicitly evaluate the quantities: 

T <0\S^{h)S\t2)S%t3)...\0> . (1.16) 

(For a detailed discussion of this approach and some third order calculations, see, for example, Q.) Since the non- 
interacting part of the Hamiltonian is independent of S, these products are actually time- independent, up to some 
minus signs arising from the time-ordering. For instance, for the S* = 1/2 case: 

r < o|5"(ti)5^(i2)|o >= 9ih - t2)S'''sy + e{t2 - h)sys^ = sgn(ii - t2)iS'. (i.i?) 

Here sgn(i) is the sign function, ±1 for t positive or negative respectively. Using the spin commutation relations 
and — S{S + 1) it is possible to explicitly evaluate the expectation values of any of the spin products occuring 
in perturbation theory. We are then left with standard fermion propogators. These are simplified by the fact that 
all fermion fields occuring in the Kondo interaction arc at r = 0; we suppress the spatial labels in what follows. For 
instance, the second order diagram is: 

-^Jdt dt'TiS^{t)S\t')) ■ T[V;t(Oy VW^ni')y^(t')], (1-18) 
which can be reduced, using Wick's theorem to: 

-J dt d^V^ 2^ ■ ^ ^^^^^ ^ < 0|V'(OV'^(i')|0 > ■ (1-19) 
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The free fermion propogator is simply: 



G{t) = (1.20) 



This gives a correction to the effective couphng constant: 

SX = ^J (1.21) 

Integrating symmetrically, Eq. (|1.19[) would give zero if the factor sgn (i), coming from the impurity spin Green's 
function, were absent. This corresponds to a cancellation between particle and hole contributions. This is as it should 
be. If we have a simple non-magnetic scatterer, with no dynamical degrees of freedom, there is no renormalization. 
The Kondo problem arises entirely from the essentially quantum-mechanical nature of the impurity spin. Including 
the sgn {t) factor, the integral in Eq. (|1.19p is infrared and ultra-violet log-divergent. In an RG transformation, we 
only integrate over a restricted range of wave-vectors, integrating out modes with D' < \k\ < D. We then obtain: 

6\ = X^h-iiD/D'). (1.22) 

The corresponding /3-function, to this order is then: 

Solving for the effective coupling at scale D, in terms of its bare value Aq at scale Dq we obtain: 

If the bare coupling is ferromagnetic, Xq < 0, then X{D) is well-behaved, getting smaller in magnitude at lower energy 
scales. However, if it is antiferromagnetic, X{D) continues to increase as we reduce the energy scale until it gets so 
large that lowest order perturbation theory for the /9-function breaks down. We may estimate the energy scale where 
this happens as: 

Tk^ Doexp{-l/Xo). (1.25) 

The scale Dq, which plays the role of the ultra-violet cut-off, is of order the band width or Fermi energy and Dq = vp/A 
where A is the cut-off in momentum units. 

After many years of research by many theorists, a very simple picture immerged for the low energy behaviour 
of the Kondo model, due in large part to the contributions of PW Anderson, ^lO] K Wilson, 11] P Nozieres[12] and 
collaborators. We may think of A as renormalizing to oo. What is perhaps surprising is that the infinite coupling 
limit is actually very simple. To see this, it is very convenient to consider a lattice model, 

oo ^ 

H = -tJ2ii'^i'^+l + ^1+1 V^O + M^^o ■ S. (1.26) 

i=Q 

The strong coupling limit corresponds to J ^ t. It is quite easy to solve this limit exactly. One electron sits at site 
and forms a spin singlet with the impurity, which I assume to have S ~ 1/2 for now. | fri> — | 4J-t>- (Here the double 
arrow refers to the impurity spin and the single arrow to the spin of the electron at site zero.) The other electrons 
can do anything they like, as long as they don't go to site 0. Thus, we say the impurity spin is "screened" , or more 
accurately has formed a spin singlet. To understand the low energy effective Hamiltonian, we are more interested in 
what the other electrons are doing, on the other sites. If we now consider a small but finite t/J, the other electrons 
will form the usual free fermion ground state, filling a Fermi sea, but with a modified boundary condition that they 
cannot enter site 0. It is as if there were an infinite repulsion at site 0. The single particle wave- functions are changed 
from sinfc(j -|- 1) to sinfcj. In the particle-hole (PH) symmetric case of a half-filled band, kp = tt/2, the phase shift 
at the Fermi surface is 7r/2. In this one-dimensional case, we take the continuum limit by writing: 

« e^'^^^Rij) + e-^'^-^Mj)- (1-27) 
For A = 0, in the PH symmetric case, the open boundary condition for the lattice model corresponds to 



^l(O) = Vfl(O) 



(1.28) 



5 



in the continimm model, just as in D=3. On the other hand, the strong coupling BC is: 



(1.29) 



The strong coupling fixed point is the same as the weak coupling fixed point except for a change in boundary conditions 
(and the removal of the impurity). We describe the strong coupling fixed point by the conformally invariant BC of 



This simple example illustrates the main ideas of the BCFT approach to quantum impurity problems. In general, 
we consider systems whose long-distance, low energy behaviour, in the absence of any impurities, is described by a 
CFT. Examples include non-interacting fermions in any dimension, or interacting fermions (Luttinger liquids) in D=l. 
We then add some interactions, involving impurity degrees of freedom, localized near r — Q. Despite the complicated, 
interacting nature of the boundary in the microscopic model, the low energy long distance physics is always described 
by a conformally invariant BC. The impurity degrees of freedom always either get screened or decouple, or some 
combination of both. Why should this be true in general? Some insight can be gained by considering the behaviour of 
arbitrary Green's functions at space- (imaginary) time points zi — ti -l-iri, Z2 = T2 + ir2, ■ ■ ■ Very close to the boundary 
we expect non- universal behaviour. If all points ri are very far from the boundary, ^ \zj — Zk\, then we expect 
to recover the bulk behaviour, unaffected by the boundary interactions. This behaviour is conformally invariant. 
However, if the time-separations of some of the points are larger than or of order of the distances from the boundary, 
which are themselves large compared to microscopic scales, then the boundary still affects the Green's functions. We 
expect it to do so in a conformally invariant way. We have a sort of conformally invariant termination of the bulk 
conformal behaviour, which is influenced by the universality class of the boundary, encoded in a conformally invariant 
boundary condition. Note that the RG flow being discussed here is entirely restricted to boundary interactions. The 
bulk terms in the effective Hamiltonian do not renormalize, in our description; they sit at a bulk critical point. We 
do not expect finite range interactions, localized near r = to produce any renormalization of the bulk behaviour. 
All of the RG flows, which play such an important role in these lectures, are boundary RG flows. 

It actually turns out to be extremely important to go slightly beyond merely identifying the low energy fixed point, 
and to consider in more detail how it is approached, as the energy is lowered. As is usual in RG analyses, this is 
controlled by the leading irrelevant operator (LIO) at this fixed point. This is a boundary operator, defined in the 
theory with the conformally invariant boundary condition (CIBC) characteristic of the fixed point. It is important to 
realize that, in general, the set of boundary operators which exist depends on the CIBC. 

In the case at hand, the simplest version of the Kondo model, the boundary operators are simply constructed out of 
the fermion fields, which now obey the BC of Eq. (|1.29p . It is crucial to realize that the impurity spin operator cannot 
appear in the low energy effective Hamiltonian because it is screened. Thus the LIO is constructed from fermion fields 
only. It must be SU(2) invariant. In general, the operator V'^"(0)V'q(0) might appear. This has dimension 1 and is 
thus marginal. Note that 1 is the marginal dimension for boundary operators in a (1-1-1) dimensional CFT since these 
terms in the action are integrated over time only, not space. If we restrict outselves to the PH symmetric case, then 
this operator cannot appear since it is odd under the PH transformation, ^ ip'^ . Thus we must turn to 4- fermion 
operators, of dimension 2, which are irrelevant. There are 2 operators allowed by SU(2) symmetry, J(0)^ and J(0)^, 
where the charge and spin currents are defined as: 



Here I am suppressing L, R indices. I work in the purely left-moving formalisim so all operators are left-movers. I 
will argue below that, since the Kondo interaction involves the spin degrees of freedom, the term in the effective 
Hamiltonian has a much larger coefficient that does the term. More precisely, we expect that coefficient of the 
term to be of order 1/Tk- By power counting, it must have a coefficient with dimensions of inverse energy. 1/Tx 
is the largest possible coefficient (corresponding to the lowest characteristic energy scale) that could occur and there 
is no reason why it should not occur in general. Another way of looking at this is that the low energy effective 
Hamiltonian has a reduced cut-off of order Tk (or Tk/vf in wave- vector units). The coefficient of is of order the 
inverse cut-off. By contrast, I shall argue below that the coeffient of ,P is of order 1/Dq <C 1/Tk- Thus, this term can 
be ignored. The precise value of the coefficient of is not known; but neither have I yet given a precise definition 
of Tk- Unfortunately, there are a large number of definitions of characteristic energy scales in use, referred to as Tn 
among other things. One possibility is to fix a definition of Tk from the coupling constant of the LIO. 



Eq. ^;2B- 



J = V^' 



(1.30) 




(1.31) 



The factor of 1/6 is inserted for convenience. The fact that the sign is negative has physical significance and in 
principle can only be deduced from comparison to other calculations (or experiments). Note that this Hamiltonian is 
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defined with the low energy effective BC of Eq. (|1.29p . This means that the unfolding transformation used to write 
Eq. P3T|) is changed by a minus sign. i.e. Eq. (|1.14p is modified to: 

^R{r) = -^L(-r), (r > 0). (1.32) 

With this effective Hamiltonian in hand, we may calculate various physical quantities in perturbation theory in the 
LIO, i.e. perturbation theory in 1/Tk- This will be discussed in detail later but basic features follow from power 
counting. An important result is for the impurity magnetic susceptibility. The susceptibility is: 

X(T) ^ 1 < [S^r >, (1.33) 

where St is the total spin operator including both impurity and electron spin operators. The impurity susceptibility 
is defined, motivated by experiments, as the difference in susceptibilities of samples with and without the impurity. 
In practice, for a finite density of impurities, it is the term in the virial expansion of the susceptibility of first order 
in the impurity density n.^: 

X = Xo + niXimp + ■■■ (1-34) 

Ignoring the LIO, Ximp vanishes at low T. This follows because the Hamiltonian of Eq. (|1.31|) is translationally 
invariant. A simple calculation, reviewed later, shows that, to first order in the LIO: 

X^mp -^r^- (1.35) 

This is the leading low T result; corrected by a power series in T/Tk- On the other hand, the high T result at 
T ^ Tif , in the scaling limit of small Aq, is: 

X^^, (1-36) 

the result for a decoupled impurity spin. Our RG, BCFT methods only give the susceptibility in the low T and high 
T limits. More powerful machinery is needed to also calculate it throughout the crossover, when T is of order Tk- 
Such a calculate has been done accurately using the Bethe ansatz solution, [l^, [l3| giving 

X{T) - -^JiT/TK) (1.37) 

where f{x) is a universal scaling function. The asymptotic results of Eqs. ()1.35p and ()1.36|) are obtained. While the 
RG, BCFT methods are generally restricted to low energy and high energy limits (near the RG fixed points) they have 
the advantages of relative simplicity and generality (i.e. they are not restricted to integrable models). The impurity 
entropy (or equivalently the impurity specific heat) has similar behavior, with a contribution in first order in the LIO: 

Sirnp QTk' (1.38) 

Again, Simp is a universal scaling function of T/Tk, approaching In 2, the result for a decoupled impurity, at high T. It 
is also possibly to calculate the impurity contribution to the electrical resistivity due to scattering off a dilute random 
array of impurities, using the Kubo formula. In this case, there is a contribution from the fixed point Hamiltonian 
itself, [i.e. from the modified BC of Eq. p.29p ] even without including the LIO. This modified BC is equivalent to a 
tt/2 phase shift in the s-wave channel (in the PH symmetric case). We may simply take over standard formulas for 
scattering from non-magnetic impurities, which make a contribution to the resistivity expressed entirely in terms of 
the phase shift at the Fermi energy. A phase shift of n/2 gives the maximum possibly resistivity, the so-called unitary 
limit: 

Pu = , ^""xa - (1-39) 

[e.VFv) 

The correction to the unitary limit can again be calculated in perturbation theory in the LIO. In this case the leading 
correction is second order: 



p{T) ~ Pu 



[T « Tk). (1.40) 
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Again, at T ;:!> Tr- we can calculate p{T) in perturbation theory in the Kondo interaction: 

In between a scaling function of T/Tk occurs. It has so far not been possible to calculate this from the Bethe Ansatz 
but fairly accurate results have been obtained using Numerical Renormalization Group methods. For a review see 

This perturbation theory in the LIO is referred to as "Nozieres local Fermi liquid theory" (FLT). This name is 
highly appropriate due to the close parallels with Fermi liquid theory for bulk (screened) Coulomb interactions. 



II. MULTI-CHANNEL KONDO MODEL 

In this lecture, I will generalize the BCFT analysis of the simplest Kondo model to the multi-channel case, 
will give a fairly sketchy overview of this subject here; more details are given in my previous summer school lecture 
notes 0- One now imagines k identical "channels" all interacting with the same impurity spin, preserving an SU{k) 
symmetry. In fact, it is notoriously difficult to find physical systems with such SU{k) symmetry, so the model is an 
idealization. Much effort has gone into finding (or creating) systems realizing the k = 2 case, with recent success. 
Jumping immediately to the continuum limit, the analogue of Eq. (jl.lSp 



The repeated index j is summed from 1 to fc; the spin indices are not written explicitly. The same "bare" EC, Eq. 
(II. lip is used. The RG equations are only trivially modified: 

=-X' + h^ + ... (2.2) 



d\nD 2 

The factor of k in the 0{X^) term follows from the closed loop in the third diagram in Fig. ([T]). We again conclude that 
a small bare coupling gets larger under the RG. However, for general choices of k and the impurity spin magnitude, 
S, it can be readily seen that the simple strong-coupling BC of Eq. (|1.29p does not occur at the infrared fixed point. 
This follows from Fig. If the coupling flowed to infinity then we would expect k electrons, one from each channel, 
to go into a symmetric state near the origin (at the first site in the limit of a strong bare coupling). They would 
form a total spin k/2. The antiferromagnetic coupling to the impurity of spin S would lead to a ground state of 
size 15 — k/2\. It is important to distinguish 3 cases: S < k/2, overscreened, S > k/2, underscreened and S = k/2 
exactly screened. It turns out this strong coupling fixed point is stable in the underscreened and exactly screened 
cases only. In the exactly screened case this follows from the same considerations as for S* = 1/2, = 1, discussed 
in the previous lecture, which is the simplest exactly screened case. Otherwise, further considerations of this strong 
coupling fixed point are neccessary. This effective spin will itself have a Kondo coupling to the conduction electrons, 
obeying the strong coupling BC of Eq. (|1.29p . This is clearest in the lattice model discussed in the previous lecture. 
The effective spin is formed between the impurity spin and the electrons on site 0. If J ^ t, electrons from site 1 can 
make virtual transitions onto site 0, producing a high energy state. Treating these in second order gives an effective 
Kondo coupling to site 1 with 

Jeff oc J. (2.3) 

The sign of Je// is clearly crucial; if it is ferromagnetic, then it would renormalize to zero. In this case, the strong 
coupling fixed point is stable. [Note that Jeff corresponds to the original J — s- oo.] However, if it is antiferro- 
magnetic then Je// gets larger as we lower the energy scale. This invalidates the assumption that the strong coupling 
fixed point is stable. [Jeff getting larger corresponds to J getting smaller.] It is not hard to see that Je// > for the 
overscreened case and Je// < for the underscreened case. [This calculation can be done in two steps. First, consider 
the sign of the exchange interaction between the electron spin on site and on site 1. This is always antiferromagnetic, 
as in the Hubbard model. Then consider the relative orientation of the electron spin on site and Seff- These are 
parallel in the overscreened case but anti-parallel in the underscreened case.] So, in the underscreened case, we may 
apply the results of the previous lecture with minor modifications. There is an impurity spin of size Seff in the low 
energy Hamiltonian but it completely decouples at the infrared fixed point. The overscreened case is much more 
interesting. It is a "non- Fermi liquid" (NFL). 
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FIG. 1: Feynman diagrams contributing to renormalization of the Kondo coupling constant to third order. 




FIG. 2: Formation of an effective spin at strong Kondo coupling, k = 3, s = 1 and Sgg = 1/2. 

To solve this case I introduce the idea of a conformal embedding (CE). This is actually useful for various other 
BCFT problems. It is a generalization of the idea of bosonization, a powerful technique in (1+1) dimensions. 

We start by considering a left-moving single component (no channels, no spin) fermion field with Hamiltonian 
density: 

n = li^hli^.. (2.4) 

Define the current (=density) operator, 

Mt + x) =: ViV-i : {x,t) = lim[i;l{x)iPL{x + e) - (0|Vi(a;)VL(a; + e)|0)]. (2.5) 

(Henceforth we generally drop the subscripts "L" and the time argument. The double dots denote normal ordering: 
creation operators on the right.) We will reformulate the theory in terms of currents (the key to bosonization). 

Consider: 

J{x)J{x + e) = : V'^(a;)V(a;)V^(a; + e)ip{x + e) : +[: i;^ {x)^{x + e) : + : ^{x)^'' {x + e) :]G(e) + G{ef 

G{e) = (0|V(a;)V^(a; + e)|0) = (2.6) 

By Fermi statistics the 4-Fermi term vanishes as e ^ 

: ij'' {x)ij{x)ij'' {x)i^{x) :=-: tj;'' {x)i;'< {x)i;{x)i;{x) : = 0. (2.7) 
The second term becomes a derivative, 

\im[J{x)J{x + e) + = lim — —[: 'il)\x)il){x + e) : — : ip^ {x + e) tp{x) :] 
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= 2i : ib^-^ib : 
ax 

n = -^J(a;)^ + constant. (2.8) 
in 

Now consider the commutator, [J(x), J(y)]. The quartic and quadratic terms cancel. We must be careful about the 
divergent c-number part, 



dx 



1 1 



X — y — iS X — y + i5 



= 2ni^6{x-y). (2.9) 

Here S is an ultraviolet cut-off. 

Now consider the free massless boson theory with Hamiltonian density (setting vf = 1): 

We can again decompose it into the left and right-moving parts, 

(p{x, t) = ipL{x + t) + iPr{x - t) 

{dt - dx)(fiL = d-(pL = 0, d+(pR = 

H = \{d.<pr + \{d+^f = lid.^uf + \{d+VLr (2.11) 

where 

d± = dt±d^. (2.12) 
Consider the Hamiltonian density for a left-moving boson field: 

n = \{d+VLf 

[d+ipLix),d+ipLiy)] = [^ + ^',ip + ^']=2i^5{x-y) (2.13) 

Comparing to the Fermionic case, we see that: 

Jl = Vnd+ipL = Vwd+cp, (2-14) 

since the commutation relations and Hamiltonian are the same. That means the operators are the same with appro- 
priate boundary conditions. 

This equivalence becomes especially powerful when the fermions have several components. Consider the case at 
hand with 2 spin components and k channels. Clearly we can write the free fermion Hamiltonian in this case as: 

n{x) = ^ lini ^ : V'^^Vaj : (x) : V'^"Vaj : {x + e) + constant. (2.15) 

aj 

It turns out to be very useful to use simple algebraic identifies to rewrite this in terms of charge, spin and channel 
current operators: 

J{x) = : ^^"'Vai : 

J = ^tai^^^. 

jA ^ ^f-^fT^y^^., (2.16) 
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Here the T s are generators of SU{k), i.e. a set of traceless Hermitian matrices obeying the orthonormaHty condition: 

TrT^T^ = is^B (2.17) 



and hence the completeness relation: 

A 

They obey the commutation relations: 



ch cd ^ cb cd 



(2.18) 



[T^,T^] =i^/^^^T^, (2.19) 

c 

where the numbers f^^'^' are the structure constants of SU{k). In the k = 2 case, we may choose: 



a 



^ y . (2.20) 



It is now straight forward to prove the following identity: 

The cofficients of each term are chosen so that the normal ordered products : ^p^"°"tpaafp^^^'4'i3b '■ and : tp^°"^ijjabfp^^^'4'i3a ■ 
have zero coefficients. The current operators now obey the current algebras: 

[J{x),J{y)] = AmkS'ix-y) 
[J''ix),J^{y)] = 2Tn5{x--y)e'''"'J'' + T:ikS''''5'{x-y) 
[J^ix),J^iy)] = 2Tn6{x-y)f^^^J^ + 2iTrS^^S\x~y). (2.22) 

The currents of different types (charge, spin, flavour) commute with each other. Thus the Hamiltonian is a sum of 
three commuting terms, for charge, spin and flavour, each of which is quadratic in currents and each of which is fully 
characterized by the current commutation relations. Upon including the right-moving degrees of freedom, we may 
define three independent field theories, for spin, charge and flavour with the corresponding Hamiltonians. The charge 
Hamiltonian is simply a free boson, with: 

J = V2^d+(p. (2.23) 

The spin and channel Hamiltonians are Wess-Zumion-Witten (WZW) non- linear cr-models (NLaM) [l^, [2l| labelled 
SU{2)k and SU{k)2. These can be written in terms of SU{2) and SU{k) bosonic matrix fields g^{t,x) and h^j{t,x) 
respectively. (These fields are Lorentz scalars; i.e. they have zero conformal spin.) The corresponding current 
operators can be written in a form quadratic in these matrix fields. In the particular case, k = 1, the corresponding 
SU{2)i WZW is simply equivalent to a free boson. The connection between the multi-component free fermion model 
and the sum of spin, charge and channel bosonic models is an example of a conformal embedding. All conformal towers 
in the free fermion finite size spectrum (FSS), with various EC's can be written as sums of products of conformal 
towers from the 3 constitutent models. Likewise, each local operator in the free fermion model is equivalent to a 
product of charge, spin and flavour local operators in the bosonic models. [Actually, this statement needs a little 
qualification. It is literally true for free fermion operators which contain even numbers of fermion fields and have zero 
conformal spin. It is only true for the fermion fields themselves if we are allowed to define chiral components of the 
WZW matrix fields.] 

We now adopt the purely left-moving representation of the Kondo model, in Eq. (jl.lSp . This has the advantage 
that the Kondo interaction can be written in terms of the spin current operators at the origin only. Thus, remarkably, 
the Kondo interaction is entirely in the spin sector: 



iJ.^-^/ drJirY + Xm.S + .... (2.24) 



Here the . . . represents the charge and channel parts of the Hamiltonian which are non-interacting, decoupled from 
the impurity. An immediate consequence of this spin-charge-channel separated form of the Kondo Hamiltonian is that 
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the Kondo interaction only appears in the spin sector. It is then reasonable to expect that the LIO at the Kondo fixed 
involves the spin operators only, corresponding to Eq. (|1.3ip . If we took Eq. (|2.24p at face value this would appear to 
be exactly true. In fact, since Eq. (|2.24p is only a low energy effective Hamiltonian, we can generate other operators, 
in the charge (and channel) sectors during intermediate stages of the RG. However, we expect all such operators to 
have much smaller coefficients, with the scale set by Dq rather than Tk- This is the reason we ignored the irrelevant 
operator in the previous lecture. The marginal operator which could be added to the effective Hamiltonian when 
particle- hole symmetry is broken, is now seen to be purely a charge operator: = \/2'Kkd+4>. Because it is linear in 
the charge boson its effects are easy to include and it is strictly marginal. Assuming it is small, we can simply ignore 
it. It leads to a line of fixed points. 

It is interesting to observe that this Hamiltonian can be formally diagonalized, for a special value of = 2/(fc + 2) 
by redefining the spin currents: 

r{r) = r{r) + 271(5(05"^. (2.25) 

It can readily be checked that the J°-{r) obey the same commutation relations as in Eq. (|2.22p . Furthermore, for 
\ — Xc, the interacting Hamiltonian reduces to: 

^=7^Z77^^/ dr[r{r)f + ... (2.26) 



27r(fc- 

This suggests that there might be an infrared stable fixed point of the RG at an intermediate value of A corresponding 
to a BCFT. To make this idea more quantitative, we must use the full apparatus of BCFT. See J. Cardy's lecture 
notes from this Summer School for a review of this subject. [TtI 

A central idea of Cardy's BCFT is that one should represent CI BC's by boundary states.^^ These states contain 
all low energy information about a BCFT. From them one can construct the finite size spectrum (with any pair of 
CI BC's at the two ends of a finite system), OPE coefficients, boundary operator content, and all other universal 
properties. To each GET there is a set of possible boundary states (i.e. a set of possible CIBC's). In general, a 
complete classification of all (conformally invariant) boundary states is not available. However, in the case of rational 
CET's, with a finite number of conformal towers, a complete classification is available. The boundary states are in 
one-to-one correspondance with the conformal towers, i.e. with the primary operators. One can obtain the complete 
set of boundary states, from a reference state by a process of fusion with primary operators. Our strategy for finding 
the low energy fixed point of the general Kondo models (and various other quantum impurity probems) is to first 
identify the boundary state corresponding to the trivial boundary conditions of Eq. (|l.lip . We then obtain the CIBC 
corresponding to the low energy fixed point by fusion with an appropriate primary operator. The choice of primary 
operator is inspired by the mapping in Eq. (I2.25|) . 

The conformal towers of the SU{2)k WZW model are labelled by the spin of the "highest weight state" (i.e. the 
lowest energy state) . [Ill, [l^, There is one conformal tower for each spin j with: 

j-0,l/2,l,...A:/2. (2.27) 

These primary fields have zero conformal spin and left and right scaling dimensions: 

A -4^. (2.28) 

We may associate them with the spin part of the fermion operators (V'l)"- The largest possible spin we can get this 
way, from n = k, anti-symmetrized with respect to flavour, is j — k/2. The fusion rules are: 

J ^ / = \j - /I, \j ~f\ + 1, b - / I + 2, . . . , min{j + f, k-j- f}. (2.29) 

Note that this generalizes the ordinary angular momentum addition rules in a way which is consistent with the 
conformal tower structure of the theories (i.e. the fact that primaries only exist with j < k/2). 

Based on Eq. (|2.25p . we expect that the infrared stable fixed point of the fc-channel Kondo model with a spin S 
impurity corresponds to fusion with the spin j = S primary operator, whenever S < k/2. Note that the spin quantum 
numbers of the conformal towers match nicely with the over/under screening paradigm. For S > k/2 we obtain the 
infrared stable fixed point by fusion with the maximal spin primary operator of spin k/2 . Thus the boundary state 
at the infrared ( "Kondo" ) fixed point is related to the free fermion boundary state by: [o, UBl 

< jO|Kondo >=< jO|free > (2.30) 
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Here the boundary states are expanded in the Ishibashi states corresponding to the Kac-Moody conformal towers of 
spin j and \jO > labels the ground state of the spin j conformal tower. S-'^, is the modular S-matrix[27j for SU{2)i^: 

This "fusion rule hypothesis" leads immediately to various predictions about the low energy behaviour which can be 
tested against numerical simulations, Bethe Ansatz calculations and experiments. One important comparison involves 
the finite size spectrum. With the free BC's of Eq. (jl.lip the FSS can be written as a sum of direct products of 
conformal towers from spin, charge and channel sectors, {Q,j,jc)- Here Q is the charge of the highest weight state 
(measured from the charge of the ground state) and jc is a shorthand notation for the SU (k) quantum numbers of 
the highest weight state of the SU{k)2 WZW model for the channel degrees of freedom. In the important example 
fc = 2 it corresponds literally to a second set of SU (2) "pseudo-spin" quantum numbers. To obtain the spectrum at 
the infrared fixed point, one replaces the spin-j conformal tower by a set of spin conformal towers using the SU{2)k 
fusion rules of Eq. (|2.29[) with j' replaced by S, the impurity spin magnitude, in the over and exactly screened 
cases. (In the underscreened case, S should be replaced by fc/2.) Since the full spectrum of each conformal tower 
is easily constructed, the "fusion rule hypothesis" predicts an infinite number of finite size energy levels. These can 
be compared to the results of Numerical Renormalization Group (NRG) calculations. These calculations give the 
spectrum of a finite chain of length /, with the impurity spin at one end, like the tight-binding model of Eq. (jl.26p . 
These spectra reveal an interesting cross-over behaviour. For a weak Kondo coupling and a relatively short chain 
length the spectrum is essentially that of the zero Kondo coupling model: i.e. the conformal spectrum with the BC 
of Eq. (|l.lip factored with the decoupled impurity spin. However, as the chain length increases this spectrum shifts. 
The characteristic cross-over length is 

^K = VF/TK^exp[l/Xo]. (2.32) 

For longer chain lengths the FSS predicted by the our BCFT methods is observed, for the low energy part of the 
spectrum. (In principle, the smaller the bare Kondo coupling and the longer the chain length the more states in this 
conformal BCFT spectrum are observed.) In [sO], for example, the first 6 energy levels (most of which are multiply 
degenerate) were compared, for the k = 2 case, obtaining excellent agreement. 



^(2j + l)(2/ + l) 
2 + k 



(2.31) 



A. Impurity Entropy 

We define the impurity entropy as: 

Si^p(T) = }nMS{l,T) - 5o(?,r)], (2.33) 

where Sq{1, T) is the free fermion entropy, proportional to Z, in the absence of the impurity. Note that, for zero Kondo 
coupling, i^jj-^p = ln[s(s + 1)], simply reflecting the groundstate degeneracy of the free spin. In the case of exact 
screening, {k — 2s), S'inip(O) = 0- For underscreening, 

^imp(0)-M5'(s' + l)], (2.34) 

where s' = s — fc/2. What happens for overscreening? Surprisingly, we will obtain, in general, the log of a non-integer, 
implying a sort of "non-integer groundstate degeneracy" . 

To proceed, we show how to calculate <S'jjj^p(0) from the boundary state. All calculations are done in the scaling 
limit, ignoring irrelevant operators, so that S^^.p(T) is a constant, independent of T, and characterizing the particular 

boundary condition. It is important, however, that we take the limit I oo first, as specified in Eq. (|2.33p . at fixed, 
non-zero T. i.e. we are interested in the limit, l/f3 — > oo. Thus it is convenient to use the expression for the partition 
fimction,^,|ii] Zab: 

Zab^J2< >< ''Ol^ > Xaie^'^"'^^^) -> e^'"/*''^ < yl|00 >< 00|B > . (2.35) 

a 

Here \aO > labels the groundstate in the conformal tower of the operator Oa and Xa is the corresponding character, 
c is the conformal anomaly. Thus the free energy is: 

Fab = -■KcT'^l/Q - Tin < yl|00 >< 00|B > . (2.36) 
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The first term gives the specific heat: 
and the second gives the impurity entropy: 



C = TTcTl/3 (2.37) 



5in^p = hi< A|00 >< 00|B > . (2.38) 

This is a sum of contributions from the two boundaries, 

Simp — SA + Ss- (2.39) 

Thus we see that the "groundstate degeneracy" gA, associated with boundary condition A is: 

expf^impA] ^|00 >= 9A- (2.40) 

Here we have used our freedom to choose the phase of the boundary state so that gA > 0. For our original, anti- 
periodic, boundary condition, g = 1. For the Kondo problem we expect the low T impurity entropy to be given by 
the value at the infrared fixed point. Since this is obtained by fusion with the spin-s (or k/2) operator, we obtain 
fromEq. ((23T|) . 

_ sin[^(2s + l)/(2 + fc)] 
^~S^~ sin[V(2 + fc)] " ^ ' 

This formula agrees exactly with the Bethe ansatz result [24| and has various interesting properties. Recall that in 
the case of exact or underscreening (s > k/2) we must replace s by k/2 in this formula, in which case it reduces to 
1. Thus the groundstate degeneracy is 1 for exact screening. For underscreening we must multiply g by (2s' + 1) to 
account for the decoupled, partially screened impurity. Note that, in the overscreened case, where s < k/2, we have: 

1 2s + l 1 , , 

< < 1 - -, (2.42) 



2 + fc 2 + k 2 + k 

so g > 1. In the case k oo with s held fixed, g ^ 2s + 1, i.e. the entropy of the impurity spin is hardly reduced 
at all by the Kondo interaction, corresponding to the fact that the critical point occurs at weak coupling. In general, 
for underscreening: 

l<5<2s+l. (2.43) 

i.e. the free spin entropy is somewhat reduced, but not completely eliminated. Furthermore, g is not, in general, an 
integer. For instance, for fc = 2 and s = 1/2, g = \f2. Thus we may say that there is a non-integer "groundstate 
degeneracy" . Note that in all cases the groundstate degeneracy is reduced under renormalization from the zero 
Kondo coupling fixed point to the infrared stable fixed point. This is a special case of a general result: the groundstate 
degeneracy always decreases under renormalization. This is related to Zamolodchikov's c-theorem[25j which states 
that the conformal anomaly parameter, c, always decreases under renormalization. The intuitive explanation of the 
c-theorem is that, as we probe lower energy scales, degrees of freedom which appeared approximately massless start 
to exhibit a mass. This freezes out their contribution to the specific heat, the slope of which can be taken as the 
definition of c. In the case of the "g-theorem" the intuitive explanation is that, as we probe lower energy scales, 
approximately degenerate levels of impurities exhibit small splittings, reducing the degeneracy. 

A "perturbative" proof of the g-theorem was given in 28] where RG flow between two "nearby" boundary RG fixed 
points with almost the same values of g was considered. A general proof was given in [2^ . 



B. Resistivity/Conductance 

In this subsection I consider the resistivity, due to scattering from a dilute array of /c-channel Kondo impurities. psj 
and the closely related conductance through a single /c-channel impurity. This latter quantity, in the k = 2 case, 
was recently measured in quantum dot experiments, as I discuss in the next lecture. Using the Kubo formula, these 
quantities can be expressed in terms of the single-electron Green's function. Due to the (5-function nature of the 
Kondo interaction, the exact retarded Green's function (in 1, 2 or 3 dimensions) with a single impurity at r = can 
be written as: 



G(r, r'; c^) = Go(|r - r*!, c^) -f Go(r, uj)T{uj)Ga{r' , lu). 



(2.44) 
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Here Go is the non-interacting Green's function. The function T, which depends on the frequency only, not the spatial 
co-ordinates, is known as the T-matrix. Note that I am using a mixed space-frequency representation of the Green's 
function, which is invariant under time-translations, but not space-translations. The only thing which distinguishes 
the dimensionality of space is Gq. 

In the case of a dilute random array of impurities, in D=3, the Green's function, to first order in the impurity 
concentration, rii, can be written exactly as: 

Gi\r~ r-U) = (2.45) 
Go (r-r'|,cj) -E(w) 

where the self-energy is given by: 

S(w) = niT{Lj). (2.46) 
(Translational invariance is restored after averaging over impurity positions.) The single-electron life-time is given by: 

T-^Lu) = ImJ:{uj) (2.47) 
and the finite temperature resistivity, p{T), can be expressed in terms of this life-time by the standard formula: 



2e'^k f tfp 



p(T) 3m2 J (27r)3 
where np is the Fermi distribution function: 



^dnp 



nF= } — r-^ . (2.49) 

exp[ep/T] + l ^ ^ 

At low temperatures, this integral is dominated by low energies so our field theory results can be used. A similar 
calculation, reviewed in the next lecture, expresses also the conductance through a single impurity in terms of Im T. 

Thus, our task is to calculate the electron Green's function in the low energy 1-dimcnsional effective field theory. 
In the zero temperature limit we may simply evaluate it at the Kondo fixed point. At low finite temperatures we 
consider the correction from the LIO. At the fixed point, the chiral Green's functions, < ipl^{r + iT)ipL{r' + it') >, 
< ip^fiir — iT)tli air' — it') > are unaffected by the Kondo interaction. We only need to consider < tpl^ir, T)ijj]i{r' , t') >. 
By general methods of BCFT this behaves as a two-point function of left-movers with the right-mover reflected to 
the negative axis, (— r',r'): 

< 0|Vr"(r,r)^H,^(r',r')|0 >= ^i'^^^.f' ,y (2.50) 

(T — t') + i(r + r) 

Only the constant, S'(i), depends on the particular CIBC. For instance, if the BC is of free fermion type, iPr{0) = 
e*''V'L(0), then iS'(i) — e"'. In general, can be expressed in terms of the boundary state, [l^l Since the fermion field 
has spin j — 1/2, the general expression is: 

^(1) - ^i/MlA^. (2.51) 
<00|A> ^ ' 

By comparing to the free fermion BC where — 1 and obtaining the Kondo BC by fusion, using Eq. (|2.30l) . it 
follow that 5(1) is given in terms of the modular S'-matrix: 

_ S'J^S^s _ cos[7r(2g + l)/(2 + fc)] 
'''~Sl/'S°s~ cos[V(2 + fc)] ■ ^'-''^ 
The zero temperature T-matrix can be expressed directly in terms of S^iy. 

= ^[1-5(1)]. (2.53) 

In this limit, T is independent of lu and purely imaginary. This result follows from the definition, Eq. (|2.44p of the 
T-matrix upon using the free Green's function: 

Go(r, uj^) = 27ne^"- [0{^LUn)O{r) - 0(^„)^^(-r)] (2.54) 



15 



where 9{x) is the step function together with the analytic continuation to real frequency: 

e{ujn) ^ e{6 - iuj) = 1. (2.55) 

In the exactly or underscreened case where we set S ~ k/2, Eq. (j2.52|) gives = — 1, the free fermion result with a 
7r/2 phase shift. This is the unitary limit resisitivity of Eq. (|1.39p . which I now define divided by a factor of k since 
we have k parallel channels. In general, at the non-Fermi liquid fixed points, 



P(0) = Pu 



1-5(1) 



< PU- (2.56) 



To calculate the leading corrections at low temperature (or frequency) we must do perturbation theory in the LIO. 
The LIO must be a boundary operator which exists under the CIBC's characterizing the fixed point and which, 
furthermore, respects all symmetries of the Hamiltonian. The set of boundary operators (for any CIBC) is a subset 
of the set of chiral operators in the bulk theory. This follows from the "method of images" approach to BCFT which 
expresses any local operator with left and right moving factors as a bilocal product of left-movers. In the limit where 
the operator is taken to the boundary, we may use the OPE to express it in terms of local left-moving operators. 
The set of boundary operators which actually exist, for a given CIBC is in one-to-one correspondance with the set 
of conformal towers in the finite size spectrum with the corresponding boundary condition imposed at both ends of 
a finite system. This can be obtained by "double fusion" from the operator content with free fermion BC's. The 
boundary operators with free BC's all have integer dimensions and include Kac-Moody descendants of the identity 
operator, such as the current operators. Double fusion, starting with the identity operators corresponds to applying 
Eq. (|2.29p twice starting with j = and j' — S. This gives operators of spin j = 0, 1, . . . min{25', k — 2S}. While 
this only gives back the identity operator j = for the exact or underscreened case, where S' = fc/2 it always gives 
j — 1 (and generally higher integer spins) for the overscreened case. We see from the dimensions, Eq. (|2.28p . that 
the spin-1 primary is the lowest dimension one that occurs with dimension 

(2.57) 



2 + k 



None of these non-trivial primary operators can appear directly in the eff'ective Hamiltonian since they are not 
rotationally invariant, having non-zero spin. However, we may construct descendant operators of spin 0. The lowest 
dimension spin zero boundary operator for all overscreened cases is J_i • ip where ip is the spin-1 primary operator. 
This is a first descendent, with scaling dimension 1 -I- A. This is < 2, the dimension of the Fermi liquid operator, 
which can also occur. This is the LIO in the exact and underscreened cases, Eq. (|1.3ip . Thus the effective 
Hamiltonian, in the overscreened case, may be written: 

H = Ho- J-i • 0. (2.58) 

Here Hq is the WZW Hamiltonian with the appropriate BC. As usual, we assume that the dimensionful coupling 
constant multiplying the LIO has its scale set by Tk, the crossover scale determined by the weak coupling RG. 
We may take Eq. (|2.58p as our precise definition of Tk (with the operator normalized conventionally). As in the 
Fermi liquid case, many different physical quantities can be calculated in lowest order perturbation theory in the LIO 
giving various generalized "Wilson ratios" in which Tk cancels. One of the most interesting of these perturbative 
calculations is for the single-fermion Green's function, giving the T-matrix. In the Fermi liquid case, the first order 
perturbation theory in the LIO gives a correction to the T-matrix which is purely real. Only in second order do we 
get a correction to Im T, leading to the correction to the resistivity of 0(l/r|-) in Eq. (|1.40p . On the other hand, 
a detailed calculation shows that first order perturbation theory in the non- Fermi liquid LIO of Eq. (|2.58p , gives a 
correction to the T-matrix with both real and imaginary parts and hence a correction to the resistivity of the form: 



p(T) = Pu 



1-5(1) 



Tk) 



(2.59) 



Here a is a constant which was obtained explicitly from the detailed perturbative calculation, having the value 
a = 4-^7? for the 2-channel S' = 1/2 case (for which = 0). Also note that the sign of the coupling constant in 
Eq. (|2.58p is not determined a priori. If we assumed the opposite sign, the T-dependent term in the resisitivy, Eq. 
(I2.59P would switch. It is reasonable to expect this negative sign, for a weak bare coupling, since the resistivity is 
also a decreasing function of T at T ^ Tk where it can be calculated perturbatively in the Kondo coupling. An 
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assumption of monotonicity of p{T) leads to the negative sign in Eq. (|2.58p . In fact, this negative sign has recently 
been confirmed by experiments, as I will discuss in the next lecture. 

A number of other low energy properties of the non-Fermi liquid Kondo fixed points have been calculated by these 
methods, including the T-depcndence of the entropy and susceptibility and space and time dependent Green's function 
of the spin density, but I will not take the time to review them here. 

III. QUANTUM DOTS: EXPERIMENTAL REALIZATIONS OF ONE AND TWO CHANNEL KONDO 

MODELS 

In this lecture I will discuss theory and experiments on quantum dots, as experimental realizations of both single 
and two channel Kondo models. 



A. Introduction to quantum dots 

Experiments on gated semi-conductor quantum dots begin with 2 dimensional electron gases (2DEG's) in semi- 
conductor heterostructures, usually GaAs-AlGaAs. (These are the same types of semi-conductor wafers used for 
quantum Hall effect experiments.) A low areal density of electrons is trapped in an inversion layer between the two 
different bulk semiconductors. Great effort goes into making these 2DEG's very clean, with long scattering lengths. 
Because the electron density is so low compared to the inter-atomic distance, the dispersion relation near the Fermi 
energy is almost perfectly quadratic with an effective mass much lower than the free electron's. The inversion layer is 
located quite close to the upper surface of the wafer (typically around 100 nm. below it). Leads are attached to the 
edges of the 2DEG to allow conductance measurements. In addition several leads are attached to the upper surface 
of the wafer, to apply gate voltages to the 2DEG, which can vary on distances of order .1 microns. Various types 
of quantum dot structures can be built on the 2DEG using the gates. A simple example is a single quantum dot, a 
roughly circular puddle of electrons, with a diameter of around .1 ji. The quantum dot is separated from the left and 
right regions of the 2DEG by large electroscatic barriers so that there is a relatively small rate for electrons to tunnel 
from the dot to the left and right regions of the 2DEG. In simple devices, the only appreciable tunnelling path from 
electrons from left to right 2DEG regions is through the quantum dot. Because the electron transport in the 2DEG is 
essentially ballistic, the current is proportional to the voltage difference Vsd (source-drain voltage) between the leads, 
for small Vsd, rather than the electric field. The linear conductance, 

/ = GVsd (3.1) 

(and also non-linear conductance) is measured versus T and the various gate voltages. 

An even simpler device of this type does not have a quantum dot, but just a single point contact between the two 
leads. (The quantum dot devices have essentially two point contacts, from left side to dot and from dot to right side.) 
As the barrier height of the point contact is raised, so that it is nearly pinched off, it is found that the conductance, 
at sufficiently low T, has sharp plateaus and steps with the conductance on the plateaus being 2ne^//i, for integer n. 
2e^ /h is the conductance of an ideal non-interacting one-dimensional wire, with the factor of 2 arising from electron 
spin. This can be seen from a Landauer approach. Imagine left and right reservoirs at different chemical potentials, 
/ifl and — sVsd, with each reservoir emitting electrons to left and right into wires with equilibrium distributions 
characterized by different chemical potentials: 

/■°° dk 

I = -2e / —v{k)[nF{ek- ^l + eVsd) - nF{ek - ^J')] 
Jo 27r 

^ „ 9 de dnp , . 2e^ , , 

G = 2e2 / ^{eu-^Ji) = —nF{-^Ji), (3.2) 
Jo 27r dfi n 

where I have inserted a factor of previously set equal to one, in the last step. Thus, provided that /i ^ fesT, 
G = 2e^//i for an ideal one-dimensional conductor. A wider non-interacting wire would have n partially occupied 
bands and a conductance of 2ne^/h. As the point contact is progressively pinched off, it is modelled as a progressively 
narrower quantum wire with fewer channels, thus explaining the plateaus. Because the gate voltage varies gradually 
in the 2DEG, backscattering at the constriction is ignored; otherwise the conductance of a single channel would be 
2e^Tr/h where is the transmission probability. 

The tunnel barriers separating the quantum dot from the left and right 2DEG regions are modeled as single 
channel point contacts. Nonetheless, as the temperature is lowered the conductance through a quantum dot often 
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tends towards zero. This is associated with the Coulomb interactions between the electrons in the quantum dot. 
Although it may be permissable to ignore Coulomb interactions in the leads this is not permissable in the quantum 
dot itself. A simple and standard approach is to add a term to the Hamiltonian of the form (5^/(2C) where Q is the 
charge on the quantum dot and C is its capacitance. In addition, a gate voltage, V, is applied to the dot, so that the 
total dot Hamiltonian may be written: 

Hi=^{h-nof, (3.3) 

where h is the number operator for electrons on the dot and no oc V. An important dimensionless parameter is t/U 
where t is the tunnelling amplitude between leads and dot. If t/?7 <C 1 then the charge on the dot is quite well-defined 
and will generally stay close to no with virtual fluctuations into higher energy states with n = no ± 1. (An important 
exception to this is when no is a half-integer. ) It is possible to actually observe changes in the behaviour of the 
conductance as ng is varied by a single step. For close to an integer value the conductance tends to become small 
at low T . This is a consequence of the fact that, for an electron to pass through the dot it must go temporarily into a 
high energy state with n = no ± 1, an effect known as the Coulomb blockade. At the special values of the gate voltage 
where no is a half-integer, the Coulomb blockade is lifted and the conductance is larger. 



B. Single channel Kondo effect 

At still lower temperatures, a difference emerges between the case where no is close to an even or odd integer; this 
is due to the Kondo effect. When no is near an odd integer the dot must have a non-zero (half-integer) spin, generally 
1/2. At energy scales small compared to U, we may disregard charge fluctuations on the dot and consider only its 
spin degrees of freedom. Virtual processes, of second order in t, lead to a Kondo exchange interaction between the 
spin on the quantum dot and the spin of the mobile electrons on the left and right side of the 2DEG. A simplified 
and well-studied model is obtained by considering only a single energy level on the quantum dot, the one nearest the 
Fermi energy. Then there are only four states available to the quantum dot: zero or two electrons or one electron 
with spin up or down. The corresponding model is known as the Anderson Model (AM): 

H = j dkiPl^'i^kaeik) +tJ dk[^l^l"d^ + h.c] + ^{nd- no)^ (3.4) 

where 

hd. = d'^°'dc,. (3.5) 

rf^" creates an electron in the single energy level under consideration on the dot. Note that I am now treating the 
conduction electrons as one-dimensional. This is motivated by the fact that the point contacts between the 2DEG 
regions and the dot are assumed to be single-channel. However, the actual wave- functions of the electrons created 
by V'fc" are extended in two dimensions on the left and right sides of the dot. The conventional label k, doesn't 
really label a wave-vector anymore. An important assumption is being made here that there is a set of energy levels, 
near the Fermi energy, which are equally spaced and have equal hybridization amplitudes, i, with the d-level on the 
quantum dot. This is expected to be reasonable for small quantum dots with weak tunnelling amplitudes, t and 
smooth point contacts. I assume, for convenience, that these wave-functions are parity symmetric between the left 
and right 2DEG's. As the gate voltage is varied, no passes through 1. Provided that T^v <C ?7, where v is the 
(1-dimensonal) density of states, we can obtain the Kondo model as the low energy effective theory at scales small 
compared to U, with an effective Kondo coupling: 

^ ~ C/(2no-l)(3-2no)" ^^"^^ 

Thus we again expect the spin of the quantum dot to be screened at T <C TV by the conduction electrons in the leads. 

A crucial, and perhaps surprising point is how the Kondo physics affects the condiictance through the dot. This 
is perhaps best appreciated by considering a tight binding version of the model, where we replace the left and right 
leads by ID tight-binding chains: 

— 2 oo jj 

H = -t J2 (4cj-Fi+^-cO-*y](4cj-Fi+/i-c-)-*'[(cli+cl)rf-F/i.c.)]-F-(nd-no)^ (3.7) 

j = -oo 1 
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The Kondo limit gives: 



H 



Cj+i + h.c.) — t 



1 



(c]cj+i + h.c.) + J(cLi + cl)-(c-i + ci) • 5 



(3.8) 



with 



J 



2t 



n 



C/(2no - 1)(3- 2no) 



(3.9) 



Note that the only way electrons can pass from the right to the left lead is via the Kondo interaction. Thus if the Kondo 
interaction is weak, the conductance should be small. The rcnormalization of the Kondo coupling to large values at 
low energy scales implies a dramatic characteristic increase of the conductance upon lowering the temperature. In the 
particle-hole symmetric case of half-filling, it is easy to understand the low temperature limit by simply taking the 
bare Kondo coupling, to infinity J ^ t. Now the spin of the quantum dot (at site 0) forms a singlet with an electron 
in the parity-symmetric state on sites 1 and —1: 



72 



-1) 



V2 



10 > 



The parity-antisymmetric orbital 



Co = 



c-i - Ci 



^/2 ' 

is available to conduct current past the screened dot. The resulting low-energy effective Hamiltonian: 



(3.10) 



(3.11) 



E (4' 



•j+i 



h.c) 



h.c.) 



x/2^^ 



+ C2)Ca + /i.e.]. 



(3.12) 



has resonant transmission, = 1 at the Fermi energy in the particle-hole symmetric case of half-filling. This leads to 

ideal 2e^/ft, conductance from the Laudaucr formula. Thus wc expect the conductance to increase from a small value 
of order at T » Tk to the ideal value at T <C Tk- Thus, the situation is rather inverse to the case of the resistivity 
due to a dilute random array of Kondo scatterers in 3 (or lower) dimensions. For the quantum dot geometry discussed 
here lowering T leads to an increase in conductance rather than an increase in resistivity. Actually, it is easy to find 
another quantum dot model, referred to as "side-coupled" where the behaviour is like the random array case. In the 
side coupled geometry the tight-binding Hamiltonian is: 



H = -t^ (ctcj+i -|-/i.c.)-h JcJ^co' 



(3.13) 



Now there is perfect conductance at J = due to the direct hopping from sites —1 to to 1. On the other hand, in 
the strong Kondo coupling limit, an electron sits at site to form a singlet with the impurity. This completely blocks 
tranmission since an electron cannot pass through without destroying the Kondo singlet. 

At arbitrary temperatures, the conductance throught the quantum dot may be expressed exactly in terms of the 
T-matrix, T{u},T). This is precisely the same function which determines the resistivity for a dilute random array of 
Kondo scatterers. To apply the Kubo formula, it is important to carefully distinguish conduction electron states in 
the left and right leads. Thus we write the Kondo Hamiltonian in the form: 



L/R 



J /• dkdk' 



27r 



i.'^i,k + ^RA'^L,k'+i^R,k')-S. 



(3.14) 



The Kondo interaction only involves the symmetric combination of left and right leads. On the other hand, the 
current operator, which appears in the Kubo formula for the conductance, is: 



3 = -^J^Wl - Nr] = -ie[H, Nl - Nr], 



(3.15) 
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where Ni^/j^ are the number operators for electrons in the left and right leads: 
Introducing symmetric and anti-symmetric combinations: 



V2 ' 



the Kondo interaction only involves ips but the current operator is: 

d f 

^ ^ lit I *['^^''^" + '^■^■l' 
which contains a product of symmetric and anti-symmetric operators. The Kubo formula: 

1 



G = hm 



e^"*<[jW,j(0)] >, 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



then then be expressed as a product of the free Green's function for -00 and the interacting one for ips. Expressing 
the Ips Green's function in terms of the T-matrix, by Eq. (|2.44p . it is not hard to show that the conductance is given 
by: 



duf 



(T) 



-7rzyImT(e, T)]. 



For T 3> Tk a perturbative calculation of the T-matrix gives: 

3^2 



2tivT 



8 



A2 + 



(3.20) 



(3.21) 



It can be checked that the higher order terms replace the Kondo coupling, A, by its renormalized value at scale to of 
T (whichever is higher), leading to the conductance: 



G 



2e' 



h 16ln' (T/Tk)' 



{T->Tk). 



(3.22) 



On the other hand, at T, a; 0, —2ttvT — > — 2i, corresponding to the 7r/2 phase shift, leading to ideal conductance. 
By doing second order perturbation theory in the LIO, Nozieres Fermi liquid theory gives: 



2e2 



Uta- 



(3.23) 



The calculation of T at intermediate temperatures and frequencies of order Tk is a difficult problem. It goes 
beyond the scope of our RG methods which only apply near the high and low energy fixed points. It is also not 
feasible using the Bethe ansatz solution of the Kondo model. The most accurate results at present come from the 
Numerical Renormalization Group method. 



C. Two channel Kondo effect 



Quite recently, the first generally accepted experimental realization of an overscreened Kondo effect, in the two- 
channel, 5* = 1/2 case, was obtained in a quantum dot device. To realize the difficulty in obtaining a two-channel 
situation, consider the case discussed above. In a sense there are two channels in play, corresponding to the left and 
right leads. However, the problem is that only the even channel actually couples to the spin of the quantum dot. In 
Eq. p.l4p it is the left-right cross terms in the Kondo interaction that destroy the two-channel behaviour. If such 
terms could somehow be eliminated, we would obtain a two-channel model. On the other hand, the only thing which 
is readily measured in a quantum dot experiment is the conductance, and this is trivially zero if the left-right Kondo 
couplings vanish. 
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A solution to this problem, proposed by Oreg and Goldhaber-Gordon[29'|, involves a combination of a small dot, in 
the Kondo regime, and a "large dot" , i.e. another, larger puddle of conduction electrons with only weak tunnelling 
from it to the rest of the system, as shown in figure ^ . This was then realized experimentally by Goldhaber-Gordon's 
group. [3^ The key feature is to adjust the size of this large dot to be not too large and not too small. It should be 
chosen to be large enough that the finite size level spacing is negligibly small compared to the other relevant energy 
scales: Tk and T which may be in the mili-Kclvin to kclvin range. On the other hand, the charging energy of the dot, 
effectively the U parameter discussed above, must be relatively large compared to these other scales. In this case, 
the charge degrees of freedom of the large dot are frozen out at low energy scales. An appropriate Kondo type model 
could be written as: 

3 3 ^ 

^ - E / dk4,i^,Mk) + |- E r.r, / ^^l^^,,.' -S+^in,- nor. (3.24) 

Here j — 1 corresponds to the left lead, j = 2 corresponds to the right lead and j = 3 corresponds to the large dot. 
Ug/i is the charging energy for the small/large dot respectivly, are the corresponding tunelling amplitudes onto the 
small dot being the total number of electrons on the large dot. no is the lowest energy electron number for the 
large dot, which is now a rather large number. (It is actually unimportant here whether no is an integer or half-integer 
because the Kondo temperature for the large dot is assumed to be negligibly small.) If Ui is sufficiently large, the 1-3 
and 2-3 cross terms in the Kondo interaction can be ignored, since they takes the large dot from a low energy state 
with na = no to a high energy state with ng, = no ± 1. Dropping these cross terms, assuming Fi = F2 as before, and 
replacing [ipi + ip2)/V2 by ips, as before, we obtain a two-channel Kondo model, but with different Kondo couplings 
for the two channels: 

= If: 

2F2 

J. ^ (3^25) 

Finally, by fine-tuning it is possible to make the two Kondo couplings equal, obtaining precisely the standard 
two-channel Kondo Hamiltonian. 

Using our BCFT methods it is easily seen that this type of channel anisotropy, with Ji 7^ J2, is a relevant 
perturbation. [30l| The relevant operator which now appears at the low energy fixed point is 

5H {Jl - J2)ipl. (3.26) 

Here ipc is a primary field in the channel sector. Since the associated channel WZW model is also SU (2)2 for k — 2 
we may label channel fields by their corresponding pseudo-spin. This primary field has pseudo-spin 1, and scaling 
dimension 1/2. It must be checked that it occurs in the boundary operator spectrum at the Kondo fixed point. This 
follows from the fact that 1^9° "P^, the product of channel and spin j = 1 primaries occurs for free fermion BC's. This 
dimension 1 operator occurs in the non-abelian bosonization formula: 

^^ri'^n^i'^'Tj^Mp^vWl (3.27) 

(Note that both sides of this equation have the same scaling dimensions and the same symmetries.) We get the 
boundary operator spectrum at the Kondo fixed point by double fusion with the j — 1/2 primary in the spin sector. 
The first fusion operation maps the js — 1 spin primary into js — 1/2 and the second one maps js — 1/2 into js — 
(and js — 1). Therefore, the jc — 1, js — primary is in the boundary operator spectrum at the (overscreened) 
Kondo fixed point. Since this operator exists and is allowed by all symmetries once the channel SU{2) symmetry 
is broken, we expect it to be generated in the low energy effective Hamiltonian. It thus destablizes the fixed point 
since it is relevant. It is not hard to guess what stable fixed point the system flows to. Suppose Ji > J2. The stable 
fixed point corresponds to J2 flowing to zero and Ji flowing to large values. The more strongly coupled channel 1 
screens the S — 1/2 impurity while the more weakly coupled channel 2 decouples. This fixed point is characterized 
by simple phase shifts of 7r/2 for channel 1 and for channel 2. Such behaviour is consistent with the weak coupling 
RG equations: 



= -Af + 7.Ai(Af + A^) + 



d\nD ' 2 

^ - -A^ + iA,(A; + A», + .... (3.28) 
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Once Af gets larger than 2A2, these equations predict that the growth of A2 is arrested and it starts to decrease. This 
RG flow is also consistent with the g-theorem: g = (1/2) In 2 at the symmetric fixed point but g — aX the stable fixed 
point. The implications of this RG flow for the conductance in the quantum dot system is also readily deduced. (3lj| 
From Eq. ((^3^ using 5'(i) = for the A: = 2, S* = 1/2 Kondo fixed point, we see that T{u! = T = 0) has half the 
value it has in the Fermi liquid case, and therefore the conductance through the quantum dot has half the value for 
the single-channel fixed point, G(0) = e^//i. Let's call the Kondo coupling to the large dot J; and the coupling to the 
symmetric combination of left and right leads, Js- The T = conductance when Jg > Ji is 2e^/h, due to the tt/2 
phase shift in the T-niatrix for the s channel. On the other hand, if Ji > J^, the T = conductance is zero since the 
s-channel phase shift is zero. For bare couplings that are close to each other, Ai w A2 the system will flow towards 
the NFL critial point, before diverging from it at low T. i.e. there is a "quantum critical region" at finite T for Ai 
near A2. The basic scaling properties follow from the fact the relevant operators destabilizing the NFL critical point 
has dimension 1/2, together with the fact that the LIO at the NFL critical point has dimension 3/2. Right at the 
critical point, the finite T correction to the conductance, of first order in the LIO is: 

G(r)^ ^[l-(7rr/TK)i/2]. (3.29) 

(The prefactor was also determined exactly here, but this is only useful if some independent measure of Tk can be 
made experimentally.) We define Tc as the crossover scale, at which the RG flow starts to deviate from the NFL 
critical point. This defines the energy scale occuring in the relevant perturbation in the effective Hamiltonian: 

H = Hnfl ± Ti/Vc(0) - -4j^J-i ■ 0s, (3.30) 

where I have included both the LIO (last term) and the relevant operator, which is present when Ai ^ A2 with a 
sign for the coupling constant cx Ai — A2. For almost equal bare couplings, this will be <C Tk, the scale at which the 
renormalized couplings become large. At Tc ^ T <C Tk, we can calculate the correction to the NFL conductance to 
first order in ip'^, giving: 

G(T) ~ -^[1 + constant • sgn(A)(Tc/T)i/2]^ (3.31) 

Here: 

A = Ai - A2. (3.32) 
We may also estimate [sT] Tc in terms of A and the average bare Kondo coupling: 

-^A^+A2^ (3.33) 

The weak coupling RG equations (to second order only) are: 



d\nD 

dA 
d\nD 



= -2AA. (3.34) 



The solution to the first of these can be written: 



for D ^ Tk- The second of these RG equations then can be written: 

dA 2 



(3.36) 



dhiD \tl{D/Tk) 
Integrating this equation gives: 

A{Tk) cx Ao/Ag, (3.37) 
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where Aq and Aq are the bare couphngs. At energy scales below Tk we may write the RG equation for A(_D): 

dA 1 , 



dbiD 2 

reflecting the fact that A has scaling dimension 1/2 at the NFL fixed point. Thus: 



(3.38) 



A{D) 



D 



1/2 



for D < Tk- By definition, the crossover scale, Tc is the energy scale where A{D) becomes of order 1. Thus: 

1/2 



1 cx AiTK) 



Tk_ 

Tc 



(3.39) 



(3.40) 



Using Eq. p.37p for A{Tk)j we finally determine the crossover scale in terms of Tk and the bare parameters: 

T.^Tk^. (3.41) 

Another interesting quantity is the dependence of the conductance, at Aq — on the temperature and the source- 
drain voltage. Vsd defines another energy scale, in addition to T so we expect: 



G 



dl 



1 - ( — j F[eV,dlT) 



(3.42) 



where F{x) is some universal scaling function. A theoretical calculation of F remains an open problem. These 
theoretical predictions, in particular the occurance of the critical exponent 1/2, are in good agreement with the 
experiments of the Goldhaber-Gordon group. [3^ 




FIG. 3: Device for realizing the 2-channel Kondo effect. 



IV. QUANTUM IMPURITY PROBLEMS IN LUTTINGER LIQUIDS 

The Kondo models considered so far in these lectures all have the property that the electrons are assumed to be 
non-interacting, except with the impurity. The validity of this approximation, is based on Fermi liquid theory ideas, 
as mentioned in the first lecture. Although our model become 1 dimensional after s-wave projection, it is probably 
important that it was originally 2 or 3 dimensional, to justify ignoring these interactions, since in 1 dimensional case, 
Fermi liquid theory definitely fails. Now interactions are important leading, at low energies, to "Luttinger liquid" 
(LL) behaviour. We will now find interesting boundary RG phenomena for a potential scatterer, even without any 
dynamical degrees of freedom at the impurity. [3^ H^l The physical applications of this theory include a point contact 
in a quantum wire, or a carbon nano-tube, a constriction in a quantum Hall bar or impurities in spin chains. 

I will just give a lightening review of LL theory here, since it has been reviewed many other places (for example 
[1, 01) and is not the main subject of these lectures. We are generally interested in the case of fermions with spin, but 
no additional "channel" quantum numbers (k = 1). A typical microscopic model is the Hubbard model: 



j 
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where hj is the total number operator (summed over spin directions) on site j. (More generally, we might consider 
"ladder" models in which case we would also get several "channels" and a plethora of complicated interactions.) 
Non-abelian bosonization is again useful, leading to a separation of spin and charge degrees of freedom. But now we 
must consider the various bulk interactions. These fall into several classes: 

• 9cJlJr- An interaction term of this form, which is proportional to {d^ipY, in the Lagrangian density, can be 
eliminated by rescaling the charge boson field: (/? \fg^- Here the Luttinger parameter, has the value 
(7 = 1 in the non-interacting case. (Unfortunately there are numerous different conventions for the Luttinger 
parameter. I follows here the notation of (45j.) This leaves the Hamiltonian in non-interacting form, but the 
rescaling changes the scaling dimensions of various operators. 

• ■0|j^^/)^^'0i|V-Lj,+h.c.: This can be bosonized as a pure charge operator. Depending on the value of the Luttinger 
parameter, it can be relevant in which case it produces a gap for charge excitations. However, this "Umklapp" 
term is accompanied by oscillating factors e*^*'^^^ so it can usually be ignored unless hp = 7r/2, corresponding to 
half-filling. It produces a charge gap in the repulsive Hubbard model at half-filling. The low energy Hamiltonian 
then involves the spin degrees of freedom only. In particular, it may correspond to the S'C/(2)i WZW model. 
In the large IJ limit of the Hubbard model, we obtain the S=l/2 Heisenberg model, with antiferromagnetic 
coupling J oc i^/C/, as a low energy E ^ U) lattice model. The low energy Hamiltonian for the Heisenberg 
model is again the SU{2)i WZW model. 

• Marginal terms of non-zero conformal spin. The only important affect of these is assumed to be to change the 
velocities of spin and charge degrees of freedom, making them, in general, different. 

• — ((7s/27r) Ji • Jfl: gs has a quadratic /3-function at weak coupling; it flows to zero logarithimically if it is initially 
positive, as occurs for the repulsive, U > Hubbard model. It is often simply ignored, but, in fact, it leads to 
important logarithmic corrections to all quantities. 

• Spin anisotropic interactions of zero conformal spin: Often SU(2) spin symmetry is a good approximation in 
materials but it is generally broken to some extent, due to spin-orbit couplings. If a U{1) spin symmetry is 
preserved then, depending on parameters, the spin degrees of freedom can remain gapless. It is then usually 
convenient to use ordinary abelian bosonization. The spin boson then also also gets rescaled ips — > gsPs where 
(7s = 1 in the isotropic case. This leads to further changes in scaling dimensions of various operators. 

• Various higher dimensional operators of non-zero conformal spin: Some of these have very interesting and 
non-trivial effects and are the subject of current research. However, these effects generally go away at low 
energies. 

Let us begin with an interacting spinless fermion model with impurity scattering at the origin only, corresponding to 
a point contact in a quantum wire. A corresponding lattice model would be, for example: 

— 1 oo oo 

H=[-t ^]V',+i-iVi^i-i^V']V'j+i+/i-c.] + C/ f^,nj+i. (4.2) 

j = -oo j=l j=~oo 

The hopping term between sites and 1 has been modified from t to i'; we might expect t' <^t for a point contact 
or constriction. Upon bosonizing and rescaling the boson, the bulk terms in the action just give: 

/oo 

dxdT{d^p^f. (4.3) 

The impurity term, in terms of continuum limit fermions, 

V2^ijj « e^'=-^Vfl(i) + e-'"->fl(j), (4.4) 

is: 

H^nt ~ ^[Jl(0) + Jfl(O) -f (4(0)Vfl(0)e''=- + h.c.)]. (4.5) 

(Note that we ignore the small variation of the continuum limit fields over one lattice spacing here. Including this 
effect only leads to irrelevant operators. This continuum limit Hamiltonian is appropriate for small \t' — t|, since we 
have taken the continuum limit assuming t' — t.) Using the bosonization formulas: 

^L/R OC e^(^±^)/v^ (4.6) 
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this becomes: 

Hint = -{f - t)V2d^e{0) - constant • (t' - t) cos[%/2(6'(0) - a)], (4.7) 

for a constant a depending on kp. While the first term is always exactly marginal, the second term, which arises 
from "backscattering" (L ^ R) has dimension 

X = g. (4.8) 

It is marginal for free fermions, where g = 1 but is relevant for g < 1, corresponding to repulsive interactions, U > 0. 
It is convenient to go a basis of even and odd channels. 

Q,l,{x) ^ ^ ' (4.9) 

The 6'e/o fields obey Neumann (N) and Dirichlet (D) BC's respectively: 

d^M - 

0o(O) - 0. (4.10) 
Then the action separates into even and odd parts, S — Se + So with: 

1 



dr / dx{d^,6eY -VbCosiOeiQ) ~ a) 







So = ^ I dT I dx{d^9oy -VfdJoiO), (4.11) 



Airg 







where Vj , the forward and backward scattering amplitudes, are both cx t' — t. The interaction term can be eliminated 
from So by the transformation: 

Ooix) ^ 9o{x) - 27TVfg ■ sgn(x). (4.12) 

On the other hand, Se is the well-known boundary sine-Gordon model which is not so easily solved. It is actually 
integrable|35l] and a great deal is known about it, but here I will just discuss simple RG results. For 5 < 1, when 
backscattering is relevant, it is natural to assume that Vb renormalizes to infinity thus changing the N boundary 
condition on 6e to D, de{0) = ct. This has the effect of severing all communication between left and right sides of 
the system, corresponding to a cut chain. If the forward scattering, Vf = then we have independent D boundary 
conditions on left and right sides: 

e{0'^)=a/V2. (4.13) 

For non-zero Vf, left and right side are still severed but the D BC's are modified to: 

e{0^) ^ a/V2TV2nVfg. (4.14) 

The simple D EC of Eq. (|4.13p or (|4.14p correspond, in the original fermion language to: 

^L(0±)cx7Afl(0±). (4.15) 

The right moving excitations on the x < axis are reflected at the origin picking up a phase shift which depends on 
Vf, and likewise for the left moving excitations on x > 0. 

Of course, we have made a big assumption here that Vf, renormalizes to 00 giving us this simple D EC. It is important 
to at least check the self-consistency of the assumption. This can be done by checking the stability of the D fixed 
point. Thus we consider the Hamiltonian of Eq. (j4.2p with t' <^ t. To take the continuum limit, we must carefully 
take into account the boundary conditions when t' = 0. Consider the chain from = 1 to cx3 with open boundary 
conditions (OEC). This model is equivalent to one where a hopping term, of strength t, to site is included but then 
a EC V'o = is imposed. From Eq. (|4.4p we see that this corresponds to '0l(O) = ~V'i?(0)- Using the bosonization 
formulas of Eq. (|4.6p we see that this corresponds to a D EC, 9{0) = constant, as we would expect from the previous 
discussion. A crucial point is that imposing a D EC changes the scaling dimension of the fermion fields at the origin. 
Setting 9{0) — constant, Eq. (|4.6p reduces to: 



(4.16) 
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The dimension of this operator is itself affected by the D BC. Decomposing ip{t,x) and 9 into left and right moving 
parts: 

f{t,x) = -^[fLit + x) + ipB.{t ~ x)] 

= V9[^L ~ ^r] (4.17) 

we see that the D BC implies: 

ipB.{0) = ¥'l(0) + constant. (4.18) 
Thus to evaluate correlation functions of (p{0) with D BC we can use: 

2 

(p(0) -> ifLiO) + constant. (4.19) 

V9 

The bulk correlation functions of exponentials of ip decay as: 

^ giai/j(t,2;)g-iaip(0,0) — g^ai^i, (t+x)/y^g-iai^i, (0,0)/v^ ^ , ^ giai/3H(*-a;)/V5g-M¥'-R.(0,0)/v^ ^ ^ 



{x + tYy^a^x-tyy^B 

(4.20) 

On the other hand, at a boundary with a D BC, 

^ gMi^(t,0)g-iai^(0,0) — <^ g2ia¥'-L(t)/%/9g-2ia(^i,(0,Q)/y^ ^ M 21) 

i(2a)V2s' ^ ■ ' 

The RG scaling dimension of the operator e"^"^ doubles at a boundary with a D BC to A = a^/g. Thus the fermion 
field at a boundary with D BC, Eq. (|4.16p . has a scaling dimension 1/25. The weak tunnelling amplitude t' in Eq. 
(|4.2p couples together two independent fermion fields from left and right sides, both obeying D BC's. Therefore the 
scaling dimension of this operator is obtained by adding the dimension of each independent fermion fields, and has 
the value 1/5. This is relevant when g > 1, the case where the weak backscattering is irrelevant, and is irrelevant 
for 5 < 1 the case where the weak backscattering is relevant. Thus our bold conjecture that the backscattering, Vf, 
renormalizes to 00 for 5 < 1 has passed an important consistency test. The infinite back-scattering, cut chain, D BC 
fixed point is indeed stable for g < 1. On the other hand, and perhaps even more remarkably, it seems reasonable to 
hypothesize that even a weak tunnelling t' between two semi-infinite chains flows to the N BC at low energies. This 
is a type of "healing" phenomena: Translational invariance is restored in the low energy, long distance limit. 

The conductance is clearly zero at the D fixed point. At the N fixed point we may calculate it using a Kubo formula. 
One approach is to apply an AC electric field to a finite region, i, in the vicinity of the point contact: 

G=lim- / dre"^^ dxT < J{y,T)J{x,0) >, (4.22) 

(independent of x). Here the current operator is J = —idrO. Using: 

< e{x, r)0(0, 0) >= -| ln(r2 + x'), (4.23) 
for the infinite length system, it is straightforward to obtain 

G = .9^. (4.24) 

(While this is the conductance predicted by the Kubo formula, it is ap par ently not neccessarily what is measured 
experimentally. [3^ For various theoretical discussions of this point see |37|.) Low temperature corrections to this 
conductance can be obtained by doing perturbation theory in the LIO, as usual. At the cut chain, D, fixed point, for 
g < 1, the LIO is the tunnelling term, cx t' in Eq. (|4.2p . This renormalizes as: 

t'{T)f^t'o{T/Tof'^-\ (4.25) 

Since the conductance is second order in t' we predict: 

G(T) cxt'oT2(i/9-i)^ ((7<l,T<To). (4.26) 
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Here Tq is the lowest characteristic energy scale in the problem. If the bare t'o is small then Tq will be of order 
the band width, t. However, if the bare model has only a weak back-scattering Vf, then Tq can be much smaller, 
corresponding to the energy scale where the system crosses over between N and D fixed points, analgous to the Kondo 
temperature. (Note that there is no Kondo impurity spin in this model, however.) On the other hand, near the N 
fixed point, where the backscattering is weak we may do perturbation theory in the renormalized back scattering; 
again the contribution to G is second order. Now, for g > I, 

Vb{T) = VboiT/nr-^ (4.27) 

and hence: 

G-geV/icx T42y2(g-i)^ (g>l,T<ro). (4.28) 

Again Tq is the lowest characteristic energy scale; it is a small cross over scale if the microscopic model has only a 
small tunnelling t' . 

A beautiful application [ssl Is^ of this quantum impurity model is to tunnelling through a constriction in a quantum 
Hall bar, as illustrated in fig. ([4]). Consider a 2DEG in a strong magnetic field at the fractional quantum Hall effect 
plateau of filling factor v — 1/3. Due to the bulk excitation gap in the Laughlin ground state there is no current 
flowing in the bulk of the sample. However, there are gapless edge states which behave as a chiral Luttinger liquid. 
Now the currents are chiral with right movers restricted to the lower edge and left movers to the upper edge in the 
figure. Nonetheless, we may apply our field theory to this system and the Luttinger parameter turns out to have 
the value g = v < 1. (Furthermore, the edge states are believed to be spin-polarized, making the spinless model 
discussed here appropriate. ) Right and left movers interact at the constriction, with a finite probability of back 
scattering, which in this case takes quasi-particles between upper and lower edges. The entire cross over function for 
the conductance can be calculated, cither by quantum Monte Carlo (together with a delicate analytic continuation to 
zero frequency) or using the integrability of the model. The results agree fairly well with experiments. 




FIG. 4: Quantum Hall bar with a constriction. Edge currents circulate clockwise and can tunnel from upper to lower edge at 
constriction. 



A number of other interesting quantum impurity problems have been studied in Luttinger liquids. These include 
the generalization of the model discussed above to include electron spin. [33| Four simple fixed points are found which 
are obvious generalizations of the two discussed in the spinless case. We may now have D or N BC for both charge and 
spin bosons, corresponding to perfect reflection/transmission for charge/spin. Interestingly, additional fixed points 
occur, for certain ranges of the charge and spin Luttinger parameters, which have charge and spin conductances 
which are universal non-trivial numbers. A general solution for these fixed points remains an open problem. A 
simpler situation occurs in spin chains, [l^ Again, I give only a telegraphic reminder of the field theory approach to 
the S=l/2 Heisenberg antiferromagnetic chain with Hamiltonian: 

H = jY^S, ■ S,+^. (4.29) 

j 

One method, is to start with the Hubbard model at 1/2-filling, where the charge excitations are gapped due to the 
Umklapp interaction. We may simply drop the charge boson from the low energy effective Hamiltonian which then 
contains only the spin boson or, equivalently, an SU {2)i WZW model. The low energy degrees of freedom of the spin 
operators occur at wave- vectors and tt: 

Sj « UMj) + + (-iy^(j), (4.30) 



27 



where the staggered component, of scahng dimension 1/2, can be written either in terms of a free boson, ip and its 
dual 6*, with 5 = 1/2 or else in terms of the primary field g'^ of the WZW model: 

/ cos((p/V2) \ 

noctrgaoc sin((p/V2) • (4-31) 
V cos(V26l) / 

The spin boson Hamiltonian contains the marginally irrelevant interaction, —{gs/2Tr)JL ■ Jr, with a bare coupling 
constant gs, of 0(1). By including a second neighbour coupling, J2, in the microscopic Hamiltonian the bare value of 
gs can be varied. At J2 = J2C ~ .2411J, a phase transition occurs with the system going into a gapped spontaneously 
dimerized phase. In the low energy effective Hamiltonian the phase transition corresponds to the bare gs passing 
through zero and becoming marginally relevant rather than marginally irrelevant. 

A semi-infinite spin chain, j > with a free BC corresponds to a D be on 9, just as for the fermionic model discussed 
above. j33|| Then the staggered spin operator at zero becomes: 

/ C0s[V2ipLm \ 

n(0)cx sin[V2^L(0)] . (4.32) 
V d^ipLiO) J 

All 3 components now have scaling dimension 1, and, it is easily seen, taking into account SU{2) symmetry that: 

n(0)(xJL(0). (4.33) 

The D BC also implies Jl(0) = Jr{0), so that both uniform and staggered spin components at a; = reduce to Jl(0). 
Now consider the effect of a Kondo type coupling between a spin chain and one additional "impurity spin" . In the 
simplest case where the impurity spin is also of size S=l/2, it makes an enormous difference exactly how it is coupled 
to the other spins. The simplest case where is where it is coupled at the end of a semi-infinite chain: 

00 

H^J'Si-S2 + jJ2 ■ ^^+1 ' (4-34) 

4=2 

with the impurity coupling J' <C J. For small J' a low energy Hamiltonian description is appropriate, and since 
S2 oc Jl(0), we obtain the contimmm limit of the Kondo model with a bare Kondo coupling A oc J'. [H, H^l Thus we 
can take over immediately all the RG results on the Kondo effect except that we must beware of logarithmic corrections 
arising from the bulk marginal coupling constant, gs, which are absent for the free fermion Kondo model. The 
correspondance with the free fermion Kondo model becomes nearly perfect when a bulk second neighbour interaction, 
J2 is added to the Hamiltonian and fine-tuned to the critical point where this bulk marginal interaction vanishes. 
Then only truly irrelevant bulk interactions (of dimension 4 or greater) distinguish the two models. The strong 
coupling fixed point of this Kondo model simply corresponds to the impurity spin being adsorbed into the chain, and 
corresponds to a renormalized J' — > J at low energies. A more interesting model involves an impurity spin coupled 
to 2 semi-infinite chains: 

-2 00 

H^jJ^Sj- Sj+i + Sj+i + J' So ■ {S-i + Si). (4.35) 

-CX3 j = l 

The continuum limit is now the 2-channel Kondo model, [33| with the left and right sides of the impurity corresponding 
to the 2 channels. Again the Kondo fixed point simply corresponds to a "healed chain" with J' renormalizing to J 
and a restoration of translational invariancc at low energies. Other possibilties involve a "side-coupled" impurity spin. 
For example we may couple the impurity spin S' to one site on a uniform chain: 

00 

H=J Sj- Sj+i + J'S' ■ Sq. (4.36) 

j = -oo 

Now the correspondance to the ordinary free fermion Kondo model fails dramatically because the boundary interaction 
oc J'S" • n{0) appears in the effective Hamiltonian where n is the staggered spin operator introduced in Eq. (j4.3ip . 
This is a strongly relevant dimension 1/2 boundary interaction. It renormalizes to infinity. It is easy to understand 
the low energy fixed point in this case by imagining an infinite bare J'. The the impurity spin forms a singlet with Sq 
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and the left and right sides of the chain are decoupled. The stability of such a fixed point is verifed by the fact that the 
spins at the end of the open chains, S±i have dimension 1 so that an induced weak "bridging" coupling Jf,ffS-i ■ Si 
has dimension 2 and is thus an irrelevant boundary interaction. We expect even a small J' to renormalize to such a 
strong coupling fixed point but in general the screening cloud will be spread over longer distances. Nonetheless, the 
left and right sides decouple at low energies and long distances. Other examples, including larger spin impurities, 
were discussed in [33l |. 

We may also couple an impurity spin to a Hubbard type model with gapless spin and charge degrees of freedom. 
The various situations closely parallel the spin chain case. In particular the cases of the impurity spin at the end of 
the chain or embedded in the middle still correspond essentially to the simple Kondo model. This follows because the 
D BC on both spin and charge bosons has the effect of reducing both uniform and staggered spin density operators 
at the boundary to Jl(0). This model, which can be applied to a quantum dot coupled to a quantum wire in a 
semi-conductor heterostructure, was analysed in detail in [4l|. 



V. QUANTUM IMPURITY ENTANGLEMENT ENTROPY 



Quantum entanglement entropy has become a popular subject in recent years because of its connection with 
black holes, quantum computing and the efficiency of the Density Matrix Renormalization Group method, and its 
generalizations, for calculating many body groundstates (on a classical computer). In this lecture I will discuss 



the intersection of this subject with quantum impurity physics. [42[ After some generalities, I will focus on the simple 
example of the single channel Kondo model, obtaining a novel perspective on the nature of the Kondo groundstate and 
the meaning of the characteristic length scale ■ In the second lecture, I discussed and defined the zero temperature 
impurity entropy, showing that it is a universal quantity, characterising the BCFT fixed point, and always decreasing 
under boundary RG flow. Quantum entanglement entropy is, in general, quite distinct from thermodynamic entropy, 
being a property of a quantum ground state and depending on an arbitrary division of a system into two different 
spatial regions. Nonetheless, as we shall see the thermodynamic impurity entropy, in the T = limit, also appears as 
a term in the entanglement entropy, in a certain limit. 

Consider first a CFT with central charge c on a semi-infinite interval, a; > 0, with a CIBC, labelled A, a.t x — 0. We 
trace out the region, cc' > a; to define the density matrix, and hence the entanglement entropy, <S'^(a;) = — trplnp, for 
the region, < x' < x. [Note that I am using the natural logarithm in my definition of entanglement entropy. Some 
authors define S using the logarithm base 2, which simply divides S by In 2.] Calabrese and Cardy (C&C) showed. [43| 
generalizing earlier results of Holzhey, Larsen and Wilczek [4^ , that this entangement entropy is given by: 

Sa{x) = (c/6) \n{x/a) + ca- (5.1) 

Here a is a non-universal constant, ca is another constant which could have been adsorbed into a redefinition of a. 
However, 5yi(x) is written this way because, by construction, the constant a is independent of the choice of CIBC, A, 
while the constant ca depends on it. C&C showed that the generalization of Sa{x) to a finite inverse temperature, 
/3, is given by a standard conformal transformation: 

Sa{x, (3) = (c/6) ln[(/3/7ra) sinh(27ra;//3)] + ca- (5.2) 

Sa{x,/3) is defined by beginning with the Gibbs density matrix for the entire system, e~^^ and then again tracing 
out the region x' > x. Now consider the high temperatures, long length limit, /? <C x: 

Sa 2ttcx/I3 + (c/6) ln(/3/27ra) + ca + 0{e-^'"'/f'). (5.3) 

The first term is the extensive term (proportional to x) in the thermodynamic entropy for the region, < x' < x. The 
reason that we recover the thermodynamic entropy when x ^ /? is because, in this limit, we may regard the region 
x' > X as an "additional reservoir" for the region < x' < x. That is, the thermal density matrix can be defined 
by integrating out degrees of freedom in a thermal reservoir, which is weakly coupled to the entire system. On the 
other hand, the region x' > x is quite strongly coupled to the region x' < x. Although this coupling is quite strong, 
it only occurs at one point, x. When x ^ /3, this coupling only weakly perturbs the density matrix for the region 
x' < X. Only low energy states, with energies of order 1/x and a neglible fraction of the higher energy states (those 
localized near x' = x) are affected by the coupling to the region x' > x. The thermal entropy for the system, with 
the boundary at x = in the limit x 3> is: 

SA,th 2ttcx/I3 + \ngA + constant, (5.4) 
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with corrections that are exponentially small in x/ j3. The only dependence on the CIBC, in this limit, is through the 
constant term, h^gA^ the impurity entropy. Thus it is natural to identify the BC dependent term in the entanglement 
entropy with the BC dependent term in the thermodynamic entropy: 

CA = ln5A- (5.5) 

This follows since, in the limit, a; 3> we don't expect the coupling to the region x' > x to affect the thermodynamic 
entropy associated with the boundary x' = 0, ca- Note that the entanglement entropy, Eq. (|5.3p . contains an 
additional large term not present in the thermal entropy. We may ascribe this term to a residual effect of the strong 
coupling to the region x' > x on the reduced density matrix. However, this extra term does not depend on the CIBC 
as we would expect in the limit x 3> /? in which the "additional reservoir" is far from the boundary. Now passing to 
the opposite limit /3 ^ oo, we obtain the remarkable result that the only term in the (zero temperature) entanglement 
entropy depending on the BC is precisely the impurity entropy, Ihqa- 

Since the impurity entropy. In 5^, is believed to be a universal quantity characterizing boundary RG fixed points, it 
follows that the boundary dependent part of the (T = 0) entanglement entropy also enjoys this property. In particular, 
we might then expect this quantity to exhibit an RG crossover as we increase x. That is, consider the entanglement 
entropy, S{x), for the type of quantum impurity model discussed in these lectures, that is described by a conformal 
field theory in the bulk (at low energies) and has more or less arbitrary boundary interactions. As we increase x, we 
might expect S{x) to approach the CFT value, Eq. (|5.1[) . with the value of ca corresponding to the corresponding 
CIBC. More interestingly, consider such a system which is flowing between an unstable and a stable CIBC, A and 
B respectively, such as a general Kondo model with weak bare couplings. Then as discussed in lectures 1 and 2, the 
impurity part of the thermodynamic entropy crosses over between two values hygA and Ings- In the general fc-channel 
Kondo case, qa = (25' + 1), the degeneracy of the decoupled impurity spin. gB is determined by the Kondo fixed 
point, having the values: 

gs = 25' + 1, (5' = 5 - 2fc, A: < 25) 
_ sin[7r(2s + l)/(2 + fc)] 



sin[7r/(2 + k)] 



[k > 25). (5.6) 



This change in lug a ought to be measureable in numerical simulations or possibly even experiments. 

To make this discussion more precise, it is convenient to define a "quantum impurity entanglement entropy" as the 
difference between entanglement entropies with and without the impurity. Such a definition parallels the definition of 
impurity thermodynamic entropy (and impurity susceptibility, impurity resistivity, etc.) used in lectures 1 and 2. To 
keep things simple, consider the Kondo model in 3 dimensions and consider the entanglement of a spherical region 
of radius r containing the impurity at its centre, with the rest of the system (which we take to be of infinite size, for 
now) . This entanglement entropy could be measured both before and after adding the impurity, the difference giving 
the impurity part. Because of the spherical symmetry it is not hard to show that the entanglement entropy reduces 
to a sum of contributions from each angular momentum channel, {I, m). For a (5-function Kondo interaction, only the 
s-wave harmonic is affected by the interaction. Therefore, the impurity entanglement entropy is determined entirely 
by the s-wave harmonic and can thus be calculated in the usual ID model. Thus we may equivalently consider the 
entanglement between a section of the chain, < r' < r, including the impurity, with the rest of the chain. It is even 
more convenient, especially for numerical simultations, to use the equivalent model, discussed in the previous lecture, 
of an impurity spin weakly coupled at the end of a Heisenberg 8=1/2 spin chain, Eq. (|4.34[) . Region A is the first r 
sites of the chain, which in general has a finite total length R. The corresponding entanglement entropy is written as 
S{ J'i^,r, R) where we set J = 1 and now refer to the impurity (Kondo) coupling as J'^. The long distance behavior 
of 5(r) is found to have both uniform and alternating parts: 

5(r) = Su{r) + {-lySAir), (5.7) 

where both Sjj and Sa are slowly varying. I will just focus here on Sjj which we expect to have the same universal 
behaviour as in other realizations of the Kondo model, including the 3D free fermion one. We define the impurity 
part of 5 precisely as: 

5™p(Jj^,r,i?) = 5c/(J^,r,i?)-5[/(l,r-l,i?-l), (r > 1). (5.8) 

Note that we are subtracting the entanglement entropy when the first spin, at site j = 1, is removed. This removal 
leaves a spin chain of length i? — 1 with all couplings equal to 1. After the removal, region A contains only r — 1 sites. 
If we start with a weak Kondo coupling, J' <g; J = 1, we might expect to see cross over between weak and strong 
coupling fixed points as r is increased past the Kondo screening cloud size S^k- Ultimately, this behaviour was found 
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but with a surprising dependence on whether R is even or odd. [Note that is a separate effect from the dependence 
on whether r is even or odd which I have aheady removed by focussing on the uniform part, defined in Eq. (|5.8p .] 

As usual, I first focus on the behaviour near the fixed points, where perturbative RG methods can be used. Strong 
coupling perturbation theory in the LIO at the Kondo fixed point turns out to be very simple. The key simplifying 
feature is that the Fermi liquid, LIO is proportional to the energy density itself, at r = 0. It follows from Eq. (|1.3ip 
and (I2.21|) (in the single channel k — 1 case) that the low energy effective Hamiltonian including the LIO at the 
Kondo fixed point can be written, in the purely left moving formalism: 



H = — dxJLixf -^Jl{0)^ dxnix)-TT^Kn{0), (5.9) 



where Ti{x) = {1 / Qn) J]^{x) is the energy density. (I have set w = 1 and used = 1/Tk-) The method of Holzhey- 
Wlczek and Calabrese-Cardy for calculating the entanglement entropy is based on the " replica trick" . That is to say, 
the partition function, Z„ is calculated on an n-sheeted Reimann surface, TZm with the sheets joined along region A, 
from r' = to r' = r. The trace of the n*'* power of the reduced density matrix can be expressed as: 

Trp(r)" = Ml (5.10) 

where Z is the partition function on the normal complex plane, C. The entanglement entropy is obtained from the 
formal analytic continuation in n: 

S = ~ \imJ-[TTp{rn (5.11) 
n^i an 

The correction to Zn{r) of first order in S,k is: 

/oo 
dr < W(r, 0) >i^„ . (5.12) 
-oo 

C&C showed related this expectation value of the energy density on TZn to a 3-point correlation function on the 
ordinary complex plane: 

^-w^^nw < H(T,0)v?„(r)(p_„(-r) >c 

<rt(r, 0)>7^„= — . (5.13) 

< V3„(r)93_„(-r) >c 

Here the primary operators ^p±n sit at the branch points ±r (in the purely left moving formulation) and have scaling 
dimension 

A„ = (c/24)[l-(l/n)2] (5.14) 
where, in this case, the central charge is c = 1. Thus Eq. (|5.13p gives: 

.'^( nw [l-(l/»)^] {2^r? 



Doing the r-integral in Eq. (|5.12p gives: 



5Z^ = -^-^n[l^{l/nn (5.16) 



Since there is no correction to Z, to first order in S^k, inserting this result in Eqs. ()5.10|) and (|5.11|) gives: 

Simp — ■ (5-17) 

Here we have used the fact that = 1, ca = at the Kondo fixed point. C&C also observed that the entanglement 
entropy for a finite total system size R, can be obtained by a conformal transformation. For a conformally invariant 
system this generalizes Eq. (|5.1|) to: 



5(r, R) = (c/6) ln[(i?/7ra) sin(7rr/i?)] + ca- 



(5.18) 
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Our result for Simp can also be extended to finite R by the same conformal transformation of the 2-point and 3-point 
functions in Eq. (j5.13p . The 3-point function now has a disconnected part, but this is cancelled by the correction to 
Z" of first order in leaving: 



_ CKTm[l - (l/n)2 



dT 



(■n/R) sinh 2iTTr/R 



sinh[7r(r + ir)/ R\ sinh[7r(r — ir)/R] 



(5.19) 



Doing the integral and taking the rephca hmit now gives: 



K 



12R 



l+TT 1- 



R 



cot 



l"!Tr\ 



(5.20) 



We emphasize that these results can only be valid at long distances where we can use FLT, i.e. r ^ £^k- They 
represent the first term in an expansion in ^a/r. 

The thermodynamic impurity entropy is known to be a universal scaling function of T/Tk- It seems reasonable 
to hypothesize that the impurity entanglement entropy (for infinite system size) is a universal scaling function of 
r/^K- At finite R, we might then expect it to be a universal scaling functions of the two variables r/^K and r/R. 
Our numerical results bear out this expectation with one perhaps surprising feature: While the scaling function is 
independent of the total size of the system at i? cx) there is a large difference between integer and half-integer total 
spin (i.e. even and odd R in the spin chain version of the Kondo model) for finite R; i.e. we must define two universal 
scaling functions Simp,e{,r/£,K, r/R) and Simp,o{r/£,K, f/R) for even and odd R. These become the same at r/i? — > oo. 

We calculated the impurity entanglement entropy numerically using the Density Matrix Renormalization Group 
(DMRG) method. In this approach a chain system is built up by adding pairs of additional sites near the centre of 
the chain and systematically truncating the Hilbert Space at a manageable size (typically around 1,000) at each step. 
The key feature of the method is the choice of which states to keep during the truncation. It has been proven that the 
optimium choice is the eigenstates of the reduced density matrix with the largest eigenvalues. Since the method, by 
construction, calculates eigenvalues of p^, it is straightforward to calculate the corresponding entanglement entropy. 
Some of our results for Simpir, R, J'j^) are shown in Figs. ([5]) and ([6]). The second figure shows that the large r result 
of Eq. (|5.20p works very well. This is rather remarkable confirmation of the universality of the quantum impurity 
entanglement entropy since we have obtained results for the microscopic model using a continuum field theory. It 
supports the idea that the impurity entanglement entropy would be given by the same universal functions for other 
realizations of the Kondo model including the standard 3D free fermion one. Note that the even and odd scaling 
functions have very different behaviour when R < For the integer spin case, Simp increases monotonically with 
decreasing r and appears to be approaching In 2 in the limit r <C R- This is what is expected from the general 
C&C result since = In 2 at the weak coupling fixed point with a decoupled impurity spin. On the other hand, for 
the half-integer spin case Simp initially increases with decreasing r/^K but eventually goes through a maximum and 
starts to decrease. The maximum occurs when r/^K « f/R] i-e. when ~ R- 

These results can be understood heuristically from a resonating valence bond picture of the ground spin in the 
Kondo spin chain model. Let's first consider the case of integer total spin, R even. Then the ground state is a spin 
singlet. It is important to realize that when J'j^ the singlet ground state becomes degenerate with a triplet state. 
This occurs since the ground state of the other R—1 sites, j = 2, 3, . . . i? has spin-1/2, as does the impurity spin and 
the two are coupled. However, by continuity in J'j^, it is the singlet state which is considered. This singlet state has 
strong entanglement of the impurity spin with the rest of the system, despite the fact that there is no term in the 
Hamiltonian coupling them together when — 0. Any singlet state, and hence the ground state for any value of J^, 
can be written as some linear combination of "valence bond states" i.e. states in which pairs of spins form a singlet. 
(We can always restrict to "non-crossing" states such that if we draw lines connecting every pair of contracted spins, 
none of these lines cross each other.) If we consider a "frozen" valence bond state, i.e. any particular basis state, then 
the entanglement entropy is simply In 2 times the number of valence bonds from region A to B. There will always 
be a valence bond connecting the impurity spin to some other spin in the system; we refer to this as the Impurity 
Valence Bond (IVB). Intuitively, if the IVB connects the impurity to a spin outside of region A then we think of 
this as resulting in an impurity entanglement of In 2; however this picture is certainly naive because the valence bond 
basis, while complete, is not orthogonal. We may think of the typical length of the IVB as being ^^j- since the spin 
screening the impurity is precisely the one forming the IVB. This picture makes it quite clear why Simp is a decreasing 
function of r and why is the characteristic scale for its variation. 

Now consider the case of half-integer spin, R odd. In this case, when = there is zero entanglement between the 
impurity spin and the rest of the system. In this case the other R—1 sites have a spin zero ground state, decoupled 
from the impurity, which is unpaired. For R odd and any the ground state always has spin 1/2. This can again be 
represented as an RVB state but each basis valence bond state has precisely one unpaired spin, which may or may not 
be the impurity. At = is is the impurity with probability 1, but consider what happens as we increase J'j^ from 
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FIG. 5: Impurity entanglement entropy for fixed rjK plotted versus r/^A' for both cases K even and R odd. FLT predictions 
for large r/^x are shown. Inset: location of the maximum, (r/5_R-)max for odd 7? plotted versus r jR. 




r 



FIG. 6: Impurity entanglement entropy for the Kondo spin chain model calculated by DRMG compared with the Fermi liquid 
calculation of Eq. (I5.20|) . 



zero, corresponding to decreasing cx exp [constant/ J^] from infinity. The probability of having an IVB increases. 
On the other hand, the typical length of the IVB when it is present is which decreases. These two effects trade 
off to give a peak in Simp when is approximately R. 

Our most important conclusion is probably that quantum impurity entanglement entropy appears to exihibit uni- 
versal cross over between boundary RG fixed points with the size, r, of region A acting as an infrared cut-off. 



VI. Y-JUNCTIONS OF QUANTUM WIRES 

Now I consider 3 semi-infinite spinless Luttinger liquid quantum wires, all with the same Luttinger parameter, 
meeting at a Y-iunction.[45| as shown in Fig. ([7]). By imposing a magnetic field near the junction we can introduce 
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a non-trivial phase, (j) into the tunneUing terms between the 3 wires. A corresponding lattice model is: 

oo 3 3 

H = J2 + ''■^■^ + Vnn,jnn+i,j] - (f /2) ^[e*^/'^i,,-^oj-i + h.c.]. (6.1) 



n=0 j=l 



(In general we can also introduce potential scattering terms at the end of each wire.) Here j = 1, 2, 3 labels the 3 
chains cylically so that we identify j = 3 with j = 0. We bosonize, initially introducing a boson <fij{x) for each wire: 



'>Pj,L/R oc exp[i{(pj ± ^j)/\/2]. 



(6.2) 



It follows immediately from the discussion of the 2- wire case in the 4*'* lecture that f is irrelevant for g < I, the case 
of repulsive interactions. Thus wc expect the decoupled wire, zero conductance fixed point to be the stable one in that 
case. On the other hand, if g > 1 the behaviour is considerably more complex and interesting. This case corresponds 
to effectively attractive interaction between electrons, as can arise from phonon exchange. 




FIG. 7: A Y-junction with voltages and currents indicated. 



It is convenient to make a basis change, analgous to the even and odd channel introduce in Lecture 4 for the 2 wire 
junction: 



V6 



(6.3) 



and similarly for the 9i's. It is also convenient to introduce 3 unit vectors, at angles 27r/3 with respect to each other, 
acting on the ($1,^2) space: 



Ki = (-l/2,V3/2) 
K2 = {-l/2,-V3/2) 
Kz = (1,0). 



(6.4) 
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We will also use 3 other unit vectors rotated by 7r/2 relative to these ones which we write as z x Ki = {—Kiy, Ki, 
The various boundary interactions are now written in this basis: 



Ti«/ = ^t^V'2L OC e-^-*g.(i/v/3)ix/?3-egZ 



r/f = ^T^^Vil OC e-H^/V^)-xKl.Hg^V^/^«o (6.5) 

et cetera. Note that we have not yet imposed any particular BC's on the fields at the origin. A more standard 
approach would probably be to start with the F = BC, 8^(0)= constant and then allow for the possibility that 
these BC's renormalize due to the tunnelling. We call the current approach the method of "Delayed Evaluation of 
Boundary Conditions". Hower, we expect the "centre of mass" field, Oq is always pinned, 8o(0) = constant. Since $o 
carries a non-zero total charge, and hence doesn't appear in any of the boundary interactions, it would not make sense 
for any other type of boundary condition to occur in the "0" sector. We may thus simply drop the factor involving Oo 
from all the boundary interactions; it makes no contribution to scaling dimensions. However, we must consider other 
possible BC's in the 1, 2 sector, since, for g > 1 the simple D BC on 9 is not stable, as mentioned above. Another 
simple possibility would be a D BC on $: $i(0) = Ci for two constants c^. To check the stability of this BC under 
the RG we must consider the LIO. The various candidates are the tunnelling and backscattering terms in Eq. (|6.5p . 
Imposing the BC $ — constant, we can evaluate the scaling dimension of the exponential factors involving 0(0) by 
the same method used in Lecture 4. The BC implies that we should replace 0(0) by 2y^^L{0)- The dimensions of 
these operators can then be read off: 

Afi^ = g/3 

Afi^ = 4g/3, (6.6) 

et cetera. T^'", for i ^ j are the LIO's. They are relevant for g < 3. Thus we see that this cannot be the stable fixed 
point for 1 < g < 3. 

Stable fixed points, for 1 < g < 3 can be identified from the form of the Tj^J^^. Iff grows large under renormalization 
it is plausible that one or the other of this set of operators could develop an expectation value. Note that if: 

< T,-,+i >^ 0, (6.7) 
this would correspond to strong tunnelling from j to j + 1. On the other hand, if: 

< T„_i >^ 0, (6.8) 

this would correspond to strong tunnelling from j to j — 1. Breaking time reversal by adding a non-zero magnetic 
flux, (p, favours one or the other of these tunnelling paths, dependng on the sign of (j). On the other hand, in the 
time-reversal invariant case cj) — Q (or vr) we do not expect such as expectation value to develop. These fixed points 
obey the BC's: 

±K^■^0) + ^/T/3{zx K^)-Q{Q) = C, (i = 1,2,3), (6.9) 

for some constants Ci. Note that since ~ 0; these are only 2 independent BC's. Introducing left and right 

moving fields: 

$ = ^($L + $fl) 

V9 

= V9i^L~^R), (6.10) 

we may write these BC's as: 

$fl,(0) -7^$L(0) + C', (6.11) 

where C is another constant vector and 7?. is a orthogonal matrix which we parameterize as: 

^^/cose -Sinn (g^l2) 
\ sm4 cos 4 / ' ^ ^ 
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with 

^ = ±2 arctan(\/3/.g). (6.13) 

We refer to these as the "chiral" fixed points, x±^ with the + or — corresponding to the sign in Eqs. (|6.9p and (|6.13p . 
Note that if ^ = we obtain the usual D BC on 9 and if ^ = tt we obtain a D BC on $. The chiral BC's are, in a 
sense, intermediate between these other two BC's. 

As usual, we check their stability by calculating the scaling dimension of the LIO. Once we have the BC's in the 
form of Eq. I|6.1ip it is straightforward to calculate the scaling dimension of any vertex operators of the general form: 

O = exp ii^d ■ $ + ) = exp[z(a - 6) • $R + i{d+h) ■ $l] oc exp{^[7^"^(a - fo) + (a + 5)] • (6.14) 



V9 



Applying this formula to the x± fixed points we find that all the operators listed in Eq. ()6.5|) have the same dimension 

3 + 5 



Ao = i|7^(a + 5) + (a-6)|2. (6.15) 



A = ^. (6.16) 



Since A > 1 for 1 < g < 3, we conclude that the chiral fixed points are stable for this intermediate range of g, only. 
Thus we hypothesize that the system renormalizes to these chiral fixed points for this range of g whenever there is a 
non-zero flux, (jj ^ 0. However, these fixed points are presumably not allowed due to time reversal symmetry when 
= and there must be some other stable fixed point to which the system renormalizes. This fixed point, which we 
referred to as "M" appears to be of a less trivial type than there "rotated D" states. An attractive possibility is that 
the M fixed point is destablized by an infinitesimal flux leading to a flow to one of the chiral fixed points. Alternatively, 
it is possible that there is a critial value of the flux, {(pd ^ neccessary to destabilize the M fixed point. In that case 
there would presumably be two more, as yet undetermined, CI BC's corresponding to these critical points. 

A. Conductance 

Once we have identified the CIBC's it is fairly straightforward to calculate the conductance using the Kubo formula. 
For a Y-junction the conductance is a tensor. If we apply voltages Vi far from the junction on lead i and let h be the 
current, flowing towards the junction, on lead i, then, for small Vi, 

3 

J.=^G.,1/,. (6.17) 

i=i 

Since there must be no current when all the Vi are equal to each other and since the total current flowing into the 
junction is always zero, it follows that 



^Gy=^Gy=0. (6.18) 



Also taking into account the symmetry of the model, we see that the most general possible form of Gij is 

— (3% - 1) + — 



G., = ^(3%-l) + -^e.,. (6.19) 



Here is the 3x3 anti-symmetric, Za-symmetric tensor with ei2 = 1. Ga, which is odd under time reversal, can 
only be non-zero when there is a non-zero flux. 

In the non- interacting case, V — in Eq. (|6.ip . g = 1, we may calculate the conductance by a simple generalization 
of the Landauer formalism. We imagine that the three leads are connected to distant reservoirs at chemical potentials 
fi — eVj . Each reservoir is assumed to emit a thermal distribution of electrons down the lead and to perfectly adsorb 
electrons heading towards it. The conductance can then be expressed in terms of the S'-matrix for the Y-junction. 
This is defined by solutions of the lattice Schroedinger equation: 

- t{<i>n+lj + *n-lj) = E^n,j, (" > 1) 

-i$i,j-(f/2)(e'^/3^o.,-i+e-'^/'$o,,+i) = Ecbo.j. (6.20) 
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The wave-functions are of the form: 

*"..An,,e-''="+Aout,,e^'", (6.21) 
for aU n > and j — 1, 2, 3, with energy eigenvalues E = — 2tcos/c. The 3x3 S-matrix is defined by: 

4ut = SA^^. (6.22) 

The most general symmetric form is: 

Sij = So {i= j) 

= S+ {t=j + l). (6.23) 

In the time- reversal invariant case 5*+ = It is then easy to see that unitarity of S implies 15*^1 < 2/3. In this 
case, for any wave- vector, k, \S±\ reaches 2/3 for some value of F of 0(1). It can also be checked that when k — (p — tt 
and r = 2, 5*+ = = 5o = and \S-\ = 1. In this case an electron incident on lead j is transmitted with unit 
probability to lead j — 1. To calculate the Landauer conductance, we observe that the total current on lead j is the 
current emitted by reservoir j, minus the reflected current plus the current transmitted from leads j ± 1: 

/dk 
— i;(fc)[(l - \S,,\')npie, - eV,) - ^ \S,.,±,\^nF{ek - eV,±^)]. (6.24) 

This gives the conductance tensor of the form of Eq. (|6.19p with: 

GsM = ^(|^+P±|^-P), (6.25) 

where S± are now evaluated at the Fermi energy. Thus the maximum possible value of G5, in the zero flux case, is 
(8/9)e^/ft-, when \S±\ — 2/3. On the other hand, for non-zero flux, when S+ = = 0, G5 — —Ga — e^//i, i.e. 
Gii = —Gi^i+i = /h but Gi^-i = 0. This imphes that if a voltage is imposed on lead 1 only, a current / = {e^ /h)Vi 
flows from lead 1 to lead 2 but zero current flows into lead 3. We refer to this as a prefectly chiral conductance tensor. 
Of course, if we reverse the sign of the flux then the chirality reverses with Vi now inducing a current {e^ /h)Vi from 
lead 1 to lead 3. 

Now consider the conductance in the interacting case, for the three fixed points that we have identified. From the 
Kubo formula, we may write the DC linear conductance tensor as: 

Gjk = lim — / dre'"" / dxT < Jj{y,T)Jk{x,Q) >, (6.26) 

where Jj = —idrdj. [At zero temperature, which I consider here, it is straightforward to take the zero frequency limit, 
in imaginary (Matsubara) formulation. At finite T it is neccessary to do an analytic continuation to real frequency 
first.] We first transform from the 4>j basis to inverting the transformation in Eq. (|6.3p : 

^3^Y.''3^^^^^■ (6.27) 
$0 makes no contribution to the conductance since Oq always obeys a D BC as discussed above. Therefore: 

_ 2 1 j'oo 

G., = lim— V v,^Vku dre^^ dxT < J^{y,T)J,{x,0) >, (6.28) 

^-0 h ^^^'^^12 J-^ Jo 

where = —idrQ^ and the result is independent of y > 0. Here 

v,^ - y^{z X Kj)^ ^ -v/273^ e^.i^;. (6.29) 

To proceed we decompose: 



J = -idAL + idr^R = -Jl + Jr- 



(6.30) 
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The general type of BC of Eq. (|6.1ip allows us to regard the <I>ij^(a;),s as the analytic continuation of the ^^^{xys to 
the negative x axis: 

$B.{x) =niL{-x) + C' (6.31) 

and thus 

Jnix) ^TZJLi-x). (6.32) 
The Green's function for Jl is unaffected by the BC: 

< JL.ir + zy)JUr' + ^x) — ^ ^. (6.33) 

2[{t-t') +i{y-x)Y 

The r integral in Eq. (|6.26p gives: 

/oo 
dre'^^T < JL^iy, t)Jl,{x, 0) >= -2ttujH{x - y)e^^y-^\ (6.34) 
-oo 

Here H{x) is the Heavyside step function, often written 9{x) but I avoid that notation here since 6{x) has another 
meaning. Thus we obtain: 

dre^"^ < J^(t, y) J,(0, x) > = -27TUj[d^,H{x - y)e"(^-") + d^,H{y - x)e^<^^-y^ 

-n^,H{y + x)e-^^y+^'^ - n,.^,H{-y ~ x)e'-'^y+% (6.35) 

Observing that H(x — y) + H{y ~ x) = 1 and H{x + y) = 1, H{~x — y) = since x and y are always positive in Eq. 
(|6.26p . we obtain: 

For the D BC on 7^ = -/ so: 

2 2 

= 2g^vj ■ Vk = 2.9^ (5,^ - 1/3), (6.37) 

corresponding to Gs — (4/3)g(e^//i), Ga = 0. We observed above that this is a stable fixed point for 5 > 3, with 
Gs > 4e^//i. This exceeds the unitary bound on the conductance in the non-interacting case. An intuitive way of 
understanding why increasing g leads to enhanced transmission is that attractive interactions can lead to pairing and 
either coherent pair tunnelling, or Andreev type processes (where an incident electron on one lead reflects as a hole 
while a pair is transmitted to a different lead) could lead to enhanced conductance. 

We can now also obtain the conductance tensor for the chiral fixed points, which are stable for 1 < g < 3. Inserting 
Eqs. dnmiand (pT^ in ((05)) we obtain: 

r e2 4.9 



/i .92 + 3 
4a2 

Ga = ±-^. (6.38) 

For g = 1 this reduces to the chiral conductance tensor discussed above in the non- interacting case with Gs ^ = 
e^/h. However, for g > 1 Ga > Gs implying that a voltage on lead 1 not only leads to all current from 1 flowing 
to 2 but some additional current also flows from 3 to 2. Intuitively, we might think that, as the electrons pass from 
lead 1 to 2 they attract some electrons from lead 3. If our hypothesis, discussed above is correct that the zero flux 
"M" fixed point is unstable, then presumably an infinitesimal flux could lead to an RG flow to these stable chiral 
flxed points. Such a device would have an interesting switching property. Even a small magnetic field could switch 
the current completely from lead 2 to lead 3, at low enough temperatures and currents. 
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VII. BOUNDARY CONDITION CHANGING OPERATORS AND THE X-RAY EDGE SINGULARITY 

There are sorae situations in condensed matter physics where we are interested in the response of a system to a 
sudden change in the Hamiltonian. A well-known example is the "X-ray edge singularity" in the adsorbtion intensity 
for X-rays in a metal, plotted versus X-ray energy. The X-ray dislodges an electron from a core level. This is assumed 
to suddenly switch on a localized impurity potential which acts on the conduction electrons. Since I have argued that 
quite generally the low energy properties of quantum impurity problems are described by CIBC's we might expect 
that the low energy response to a sudden change in impurity interactions might be equivalent to the response to a 
sudden change in CIBC's. Very fortunately, Cardy also developed a theory of BC changing operators which can be 
applied to this situation. In this lecture I will show how this theory can be applied to the usual X-ray edge problem 
and to a multi-channel Kondo version. 1471 



A. The X-Ray Edge Singularity 

When an X-ray is adsorbed by a metal it can raise an electron from a deep core level, several keV below the Fermi 
surface, up to the conduction band. Let Eq be this large energy difference between the Fermi energy and the core 
level and let w be the energy of the X-ray. At T = 0, ignoring electron-electron interactions, this transition is only 
possible for uj > Eq. Here Eq is a "renormalized" value oi Eq. I am assuming that the core level has a distinct 
energy, rather than itself being part of an energy band. This may be a reasonable approximation since core levels 
are assumed to be tightly bound to nuclei and to have very small tunnelling matrix elements to neighbouring nuclei. 
Presumably the excited electron will eventually relax back to the core level, possibly emitting phonons or electron-hole 
pairs. This is ignored in the usual treatment of X-ray edge singularities. Thus we are effectively ignoring the finite 
width of the excited electron states. Thus the X-ray adsorption intensity, /(w) will be strictly zero for co < Eq, in 
this approximation. When the core electron is excited into the conduction band, it leaves behind a core hole, which 
interacts with all the electrons in the conduction band. Note that the only interaction being considered here is the one 
between the core hole and the conduction electrons. The X-ray edge singularity, at a; = i?o, in this approximation, is 
determined by the response of the conduction electrons to the sudden appearance of the core hole potential, at the 
instant that the X-ray is adsorbed. It turns out that, for uj only slightly larger then Eq, very close to the threshold, 
I{uj) is determined only by the conduction electron states very close to ep', this fact allows us to apply low energy 
effective Hamiltonian methods. The difference between Eq and Eq arises from the energy shift of the filled Fermi sea 
due to the core hole potential. Not including the interaction with the external electromagnetic field, which I turn to 
momentarily, the Hamiltonian is simply: 

+ Eob'<b, (7.1) 

where b annihilates an electron in the core level at r = 0. I have assumed, for simplicity, that the core hole potential, 
VS^ (r) , is a spherically symmetric 5- function. These assumptions can be easily relaxed. Following the same steps as 
in Sec. I, we can reduce the problem to a one-dimensional one, with left-movers only: 

^ " ^^^irf;:^^ + ^^^Vi(0)V'L(0) + EQb^b. (7.2) 

(I have set vp = 1 and V cx V.) It is convenient to bosonize. We may introduce a left-moving boson only: 

V'icxe'^"^^ (7.3) 



H 



i/;t(r) 



2m 



i^ir) + Vbb'<d^{r)tlj^ 



H = / (dM^ - ^bb^d,q>L + EQbb^. (7.4) 

The solubility of this model hinges on the fact that b^b commutes with H. Thus the Hilbert Space breaks up into 
two parts, in which the core level is either empty or occupied. The spectrum of the Hamiltonian, in each sector of 
the Hilbert Space, is basically trivial. In the sector where the core level is occupied, = 1, we get the spectrum of 
free electrons with no impurity: 

/•OO 

Hq = / (9.0i)^ (7.5) 
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In the sector with 



we get the spectrum with a potential scatterer present: 



poo y 



(7.6) 



What makes this problem somewhat non-trivial is that, to obtain the edge singularity, we must calculate the Green's 
function of the operator which couples to the electromagnetic field associated with the X-rays: il!'^{t,r = 0)b{t). This 
operator, which excites an electron from the core level into the conduction band, mixes the 2 sectors of the Hilbert 
Space. There is also some interest in calculating the Green's function of the operator b{t) itself; this is associated 
with photo-emission processes in which the core electron is ejected from the metal by the X-ray. Again this operator 
mixes the two sectors of the Hilbert space. 

In what may have been the first paper on bosonization, in 1969, Schotte and Schotte [i^ observed that these 
Green's functions can be calculated by taking advantage of the fact that the two Hamiltonians, Hq and Hi are related 
by a canoical transformation (and a shift of the ground state energy). To see this note that we may write Hi in the 
form: 



Hi 



{dx4>L)'^ + constant, 



(7.7) 



where: 



Using the commutator: 



(t>L{x) = (I>l{x) 



V 



sgn(a;). 



(7.8) 



[dy(l>Liy),<t)Liy)] 



(7.9) 



we see that: 



Hi = U^'HoU + constant, 



(7.10) 



with the canonical transformation: 



C/ = expHl/0L(O)/V^. 



(7.11) 



(U can only be considered a unitary operator if we work in the extended Hilbert space which includes states with all 
possible BCs. U maps whole sectors on the Hilbert space, with particular BCs, into each other.) Consider the core 
electron Green's function, < 5(i)t6(0) >. Wc may write: 



(7.12) 



This is valid because, due to Fermi statistics, the core level must be vacant before acts, and occupied after it acts, 
i.e. we can replace the 66^ factor in by 1 on the right hand side and by on the left. But Hq and Hi commute 
with b so we have: 



< 0|6t(i)6(0)|0 1|6^&|1 >< 0\e'"°'e-'"''\0 > 



(7.13) 



Here |0 > is the ground state of the system, including the core level and the conduction electrons. This state can be 
written: |0 >= |1 > |0 > where |1 > is the state with the core level occupied and |0 > is the filled Fermi sea ground 
state of the conduction electrons, with no impurity potential. Using Eq. (|7.10p we see that: 



< Q\\e'"°'e-'"''\b >= e-''-^°* < b\e'"°'U^ e-'^'^'UlO >= er'^"' < 0|C/(<)^C/(0)|0 > 



(7.14) 



Using our explicit expression, Eq. (|7.1ip for U, we have reduced the calculation to one involving only a free boson 
Green's function: 



< 0|6^(<)6(0)|0 >= 



>cx 



Similarly, to get the Green's function of {t)xi;L{t,Q) I use: 



b\t)i,L{m 



(7.15) 



(7.16) 
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and thus 

<0|6t(t)V;i(i,0)6(0)V^i(0,0)|0> = e-*^°* <6|t/t(t)?7(0)^L(t,0)V'I(0,0)|6> 



Finally, we Fourier transform to get the X-ray edge singularity: 

/oo i{iAj — Eo)t Oluj P ') 

This result is conventionally written in terms of a phase shift at the Fermi surface, rather than the potential 
strength, V . The connection can be readily seen from Eq. (|7.8p and the bosonization formula: 

i^L{x) (X e'v^0.(-)e*V'.sgn(x)/2^ (^^^g) 
Since "tphix) for x < represents the outgoing field, we see that: 

^out = e'^Vin, (7.20) 

where the phase shift is: 

5^-V/2. (7.21) 

In fact, the parameter V appearing in the bosonized Hamiltonian should be regarded as a renormalized one. Its 
physical meaning is the phase shift at the Fermi surface induced by the core hole. Only for small V is it linearly 
related to the bare potential. Even if the core hole potential has a finite range, we expect the formulas for the X-ray 
edge singularity to still be correct, when expressed in terms of the phase shift at the Fermi surface, 5. More generally, 
if the core hole potential is not a 5-function, but is still spherically symmetric, a similar expression arises for the 
X-ray edge singularity with the exponent involving a sum over phase shifts at the Fermi surface, 5i, in all angular 
momentum channels, I. 

The connection with a boundary condition changing operator (BCC0)[l3, [ill is now fairly evident. If we revert to 
the formulation of the model on the semi-infinite line, r > 0, then the boundary condition is: 

V-flW^e^^ViCO). (7.22) 

The operator, U or 6, which creates the core hole potential in the Hamiltonian can be viewed as changing the boundary 
condition, by changing the phase shift 5. It is interesting to consider the relationship between the finite size spectrum 
with various EC's and the scaling dimensions of h and '0i(O)''^&. We consider the system on a line of length I with 
a fixed BC iPr{1) — — V'l(0 the far end. Equivalently, in the purely left-moving formulation, for (5 = we have 
anti-periodic BC's on a circle of circumference 21: iPl{x + 21) = —iPl{x). This corresponds to periodic BC's on the 
left-moving boson field, 

0l(x + 2O = </)L(a;) + x^g, (g = 0,±l,±2,...). (7.23) 

The mode expansion is: 

(t + x) °° 1 

0L(i, x) = \/^^r7-^(3 + V , [exp(-i7rm(t -I- x)/l)am + h.c] . (7.24) 
21 ^ — \J2ixm 

m—l ^ 



The finite size spectrum is: 



oo 



1 X • 

m—l 



(7.25) 



The universal ground state energy term, — 7r/(24/) has been included. We see that Q can be identified with the charge 
of the state (measured relative to the filled Fermi sea). There is a one to one correspondance between the states in 
the FSS with excitation energy irx/l and operators with dimension x in the free fermion theory. For example, Q = ±1 
corresponds to V'l and ipj^ respectively, of dimension x = 1/2. If we impose the "same" BC at both ends: 

^'^(O) = e^^'i^LiO) 

M) = -e^^'Ml) (7.26) 
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then this corresponds to the same anti-periodic BC on tpL{x) in the purely left-moving formulation, so the spectrum 
is unchanged. On the other hand if the phase shift is inserted at a: = only, then + 21) = —e~^'^^'ipL{x), 

corresponding to: 

cj)L{^ + 2l)-c^L{0) = V^{n-S/7r), (7.27) 
corresponding to the replacement Q ^ Q — S/tt in Eq. (|7.24p . The FSS is now modified to: 



E = / {d^c^L? 



m— 1 



In particular, the change in ground state energy due to the phase shift is 

Til fS^"^ 



(7.28) 



E^{5)~E^{Q) = --^-j . (7.29) 

Actually, adding the potential scattering also changes the ground state energy by a non-universal term, of 0(1) which 
was adsorbed into Eq in the above discussion. It is the term of 0(1//) which is universal and determines scaling 
exponents. The corresponding scaling dimension, x = ((5/7r)^/2 is precisely the scaling dimension of the BCCO 6 (or 
U). The energy of the excited state with Q = \ obeys: 

E{Q,6)~En{G) = ^-(l~-^ . (7.30) 

The corresponding scaling dimension, x = (1 — 5/tt)'^ /2, is the scaling dimension of ?/'l(0)^6. 

This is all in accord with Cardy's general theory of BCCO's. An operator which changes the BC's from A to B 
generally has a scaling dimension, x, which gives the ground state energy on a finite strip of length I with BC's A at 
one end of the strip and B at the other end, measured relative to the absolute ground state energy with the same BC's 
A, at both ends. This follows by making a conformal transformation from the semi-infinite plane to the finite strip. 
Explicity, consider acting with a primary BCCO O at time ri at the edge, x = of a semi-infinite plane. Assume O 
changes the BC's from A to B. Then at time T2 change the BC's back to A with the hermitean conjugate operator 
0\ also acting at a; = 0, as shown in Fig. ((8|). Then the Green's function on the semi-infinite plane is: 

< A\Oin)OHr,)\A >^ (7.31) 

(T2 - Tl)^^ 

where x is the scaling dimension of O. Now we make a conformal transformation from the semi-infinite plane, 
z = T + ix, {x > 0) to the finite width strip: w = u + iv, < v < I: 

z = /e^"'/'. (7.32) 

Note that the positive real axis, x = 0, r > 0, maps onto the bottom of the strip, v = 0. Choosing ti, T2 > 0, both 
points map onto the bottom of the strip. Note that, on the strip, the BC's are A at all times at the upper boundary, 
V = I but change from ^ to _B and then back to A on the lower strip, at times ui, U2- The Green's function on the 
strip is given by: 

<AA\0{u,)0\u2)\AA>= { — -) . (7.33) 

\ 21 smh[(7r(wi — U2)/ [2l)\ J 

Here \AA > denotes the ground state of the system on the strip of length / with the same BC, A, at both ends of the 
strip. We may insert a complete set of states: 

< AA\0{u,)0^ {u2)\AA ^ I < AA\0\n > |2e-s„(„.-ni)^ (7 34) 

n 

However the states \n > must all be states in the Hilbert Space with different BC's A at v — and B at v = I. The 
corresponding energies. En are the energies of the states with these different BC's measured relative to the absolute 
ground state energy with the same BC's A at both ends. Taking the limit of large U2 — ui in Eqs. (|7.33p and (|7.34[) 
we see that 

T 

where is the lowest energy state with BC's A and B produced by the primary operator O. The lowest dimension 
BCCO will simply produce the ground state with BC's A and B. This corresponds to the operator 6 in the X-ray 
edge model. On the other hand, the lowest energy state produced could be an excited state with BC's A and B as in 
the example of bip]^{0), which produces the primary excited state with Q = 1 and BC's twisted by the phase 6. 



(E^-E,), (7.35) 
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FIG. 8: BCCO's act at times ri and T2 on the semi-infinite plane. This is conformally mapped to the infinite strip with the 
BCCO's acting on the lower boundary. 



B. X-ray Edge Singularites and the Kondo Model 

This BCCO approach has various other applications that go beyond the Schotte and Schotte results. One of 
them is to the Kondo model. So far we have been ignoring electron spin. If the core level is doubly occupied in the 
ground state, then it would have spin-1/2 after one electron is ejected from it by the X-ray. An initially spin-0 ion 
would then acquire a net spin-1/2. In addition to the potential-scattering interaction with the localized core hole, 
there would also be a Kondo interaction. Thus, it is interesting to consider the effect of suddenly turning on a Kondo 
interaction. We might expect that the Kondo effect could dominate the X-ray edge exponent, at least at low enough 
temperatures T <C Tk and frequencies: to — Eq <^ Tk- Thus, we consider the Hamiltonian: 

/oo J -'13 

dr^T^^^Uc + AV'r^(0)^fc5(0) • 6"t^6/5 + E^h^^h^. (7.36) 
-oo "'^ 2 2 

Now creates a core electron with spin a and we have considered the general case of k channels of conduction 
electrons, i = z, 2, 3, . . . Again we are interested in Green's functions for the core electron operator, ha{t) and also 
the operators b"'' {t)ip]^iij{t, 0). These Green's functions should exhibit a non-trivial cross-over with frequency, or time, 
but at long times (frequencies very close to the threshold) we expect to be able to calculate them using properties of 
the Kondo fixed point. Since the operator ba creates the impurity spin, thus turning on the Kondo effect, it is again 
a BCCO. In this case, it should switch the EC from free to Kondo. Thus we expect its infrared scaling dimension to 
be given by the energy of the ground state with a free BC at one end of the finite system and a Kondo EC at the 
other. This spectrum is given by fusion with the j — 1/2 primary in the spin sector. The ground state with these 
EC's is always the j = 1/2 primary itself, of dimension 

x^(3/A)/{2 + k). (7.37) 

This follows because the free spectrum includes the charge zero, spin j = fiavour singlet, (0,0,/). Fusion with 
the spin j primary always gives 0, j, I) among other operators. This appears to have the lowest dimension of all 
fusion products. Note that this operator does not occur in the operator spectrum considered earlier at the Kondo 
fixed point. There we only considered operators produced by double fusion, corresponding to the FSS with Kondo 
EC's at both ends of the finite system. This gives the operator spectrum with a fixed, Kondo EC. Eut for the 
Kondo X-ray problem, we must consider the corresponding ECCO. In general, for a CIEC obtained by fusion with 
some operator O from a non-interacting EC, we may expect that the ECCO will be O itself. We may check this 
result by a more elementary method in the single channel case. There the Kondo fixed point is equivalent to a phase 
shift of magnitude 7r/2 for both spin up and spin down, S^^^. The energy of this state, from Eq. (j7.29p is simply 
(l/27rZ)[((5|/7r)^ -|- (Si/tt)'^] implying a dimension x ~ 1/4. This agrees with Eq. (|7.37p in the special case k = 1. We 
may also consider the dimensions of the operators V'l"^(0)&/3. This operator has Q = 1 (one extra electron added 
to the conduction band), transforms under the fundamental representation of flavour, of dimension k and has spin 
either j — or j = I depending on how we sum over the spin indices a and (3. The free spectrum always contains 
the operator corresponding to the fermion field itself, {Q — l,j = 1/2, k) (where k now denotes the /c-dimensional 
fundamental representation of SU{k)) and fusion with j = 1/2 gives {Q — — 0, k) for all k and {Q = l,j = 1, fc) 
for fc > 2. These operators have dimension: 

4fc + 2fc(2Tfc) + ^T^- = ^'-''^ 

(It can be seen than 2:1/2 = 1/2 corresponding to the free fermion operator.) Again we may check the case fc = 1 by 
more elementary arguments. We may find the unitary operators corresponding to b^. as: 

6q cx exp[2i((5f(?!)|L -I- (7.39) 
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where we have introduced separate bosons for spin up and spin down electrons. can depend on a. It is 

convenient to switch to charge and spin bosons, 

0c/. (7.40) 
By choosing {5^,6i) = {tt/2,—tt/2) for 6| and (<5|,5j,) = (— 7r/2,7r/2) for 6| we obtain: 

byi oc exp(±iV2TT(t)Ls)- (7-41) 
These have the correct quantum numbers as can be seen by comparing with the standard bosonization formula 

for ^La- 

i/jyiL oc exp{iV2TT(pLc) exp{±iV2TT(pLs)- (7.42) 

exp{±i\/2n(j)Ls) can be identified with gLyi the chiral component of the WZW model fundamental field. We then 
see that the spin singlet operator, exp(i-\/27r0ic), has x = 1/4 in agreement with Eq. (|7.38|) . On the other hand, 
the triplet operators have dimension x — 5/4. Since there is no j = 1 primary for k = 1 they contain Kac-Moody 
descendents, i.e. the spin current operator. 

VIII. CONCLUSIONS 

Apart from the examples discussed in these lectures, BCFT techniques have been applied to a number of other 
quantum impurity problems, including the following. We can consider a local cluster of impurities. At distances large 
compared to the separation between the impurities the same methods can be applied. The 2-impurity Kondo model 
exhibits a NFL fixed point which can be obtainedjisj by a conformal embedding which includes an Ising sector in 
which the fusion is performed. The 3-impurity Kondo model also exhibits a novel NFL fixed point. It was obtained 49| 
by a different conformal embedding with fusion in a parafermion CFT sector. Impurities in SU (3) spin chains 50j 
and quantum Brownian motionjsij were also solved by these techniques. They were even applied [53| to a high energy 
physics model associated with Callan and Rubukaov. This describes a super-heavy magnetic monopole interacting 
with fc-flavours of effectively massless fermions (quarks and leptons) . The monopole is actually a dyon having a set of 
electric charge states as well as a magnetic charge. When the fermions scatter off the dyon they can exchange electric 
charge. In this case fusion takes place in the charge sector. 

The assumption that essentially arbitrary impurity interactions, possibly involving localized impurity degrees of 
freedom, interacting with a gapless continuum, renormalizes at low energies to a CIBC has worked in numerous 
examples. It appears to be generally valid and will likely find many other applications in the future. 
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